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A  LASER  TEST  SET  FOR  THE  LOW-POWER  ATMOSPHERIC 
COMPENSATION  EXPERIMENT  SATELLITE 


1.  INTRODUCTION 

The  spatial  profile  of  a  light  beam  altered  by  transmission  through  atmospheric  turbulence  can  be 
restored  by  compensation.  Atmospheric  compensation  techniques  rely  on  knowledge  of  the  altered  beam 
intensity  profile.  The  Low-Power  Atmospheric  Compensation  Experiment  (LACE)  satellite  was 
developed  for  proof-of-concept  tests  for  atmospheric  compensation  of  ground-based  lasers.  The  objective 
of  LACE  is  to  measure,  in  real  time,  the  absolute  intensity  and  spatial  distribution  of  low-energy, 
atmospherically  nondamaging  ultraviolet,  visible,  and  infrared  radiation  transmitted  from  ground  sites  to 
the  satellite.  An  array  of  335  sensors  is  on  the  face  of  the  satellite  that  is  used  to  measure  the  intensity 
profile  of  the  laser  radiation  incident  on  the  satellite  in  orbit.  This  array  comprises  85  sensor  pairs  for 
pulsed  ultraviolet-visible  near-infrared  signals,  85  sensors  for  acousto  optically  modulated  visible  near- 
infrared  signals,  and  40  sensor  pairs  for  mechanically  chopped  infrared  signals.  Figure  1  schematically 
shows  the  sensor  array. 

Prior  to  launch,  two  functional  tests  were  required  for  sensor  and  electronics  performance 
verification:  sensor  integration  in  a  model  spacecraft  (integration  test)  and  operational  evaluation  in  the 
completed  spacecraft  under  vacuum  and  thermal  cycling  (vacuum  test).  Because  of  the  large  number  of 
sensors  and  the  short  satellite  operating  times  imposed  by  cooling  requirements,  it  was  imperative  that 
irradiation  of  each  sensor  occur  rapidly  and  automatically.  To  conduct  these  tests,  a  laser  test  set  (LATS) 
was  developed.  The  principal  function  of  the  test  system  was  to  irradiate  the  sensor  array  subsystem 
(SAS)  with  laser  radiation  replicating  that  to  be  employed  in  the  LACE  experiments.  This  provided 
information  to  the  LACE  ground  computer  on  laser  power  levels  at  the  individual  sensors  to  verify 
wavelength  response  and  dynamic  range  of  the  SAS.  LATS  was  not  designed  to  quantitatively  calibrate 
sensor  parameters  but  to  qualitatively  verify  sensor  performance.  Both  individual  sensor  tests  (small 
cross-section  beams  scanned  over  the  face  of  the  satellite)  and  mukisensor  flood  tests  of  the  entire  satellite 
face  were  required  for  validation.  Repeated  irradiation  of  each  individual  sensor  required  a  precise  and 
reproducible  pointing  scheme.  Flood  tests  of  the  entire  sensor  array  required  a  removable  Icns/mirror 
combination  to  produce  an  expanded  beam.  The  LATS  was  used  in  the  integration  test,  and  then  it  was 
modified  for  use  in  the  vacuum  tests. 

2,  LATS 

Figure  2  is  a  schematic  drawing  of  the  LATS  configuration.  The  elements  of  LATS  included  a  set 
of  lasers,  a  transfer  mirror,  a  goniometer-mounted  steering  mirror  and  stepping-motor  controller,  a  power 
meter,  a  computer,  and  a  network  of  bidirectional  data  buses.  The  laser  beam  was  expanded  and 
collimated  to  a  2-in.  diameter  to  provide  a  uniform  beam  profile.  To  provide  an  intensity  reference,  a 
beamsplitter  reflected  a  portion  of  the  light  onto  a  power  meter.  The  computer  recorded  the  reference 
power  by  using  an  analog-to-digital  (A/D)  converter.  The  transmitted  portion  of  the  laser  beam  was 
reflected  from  a  transfer  mirror  (TM)  located  in  front  of  the  plane  of  the  sensor  panel  onto  the 
goniometer-mounted  steering  mirror  (SM).  The  computer  then  aimed  the  steering  mirror  to  direct  the 
laser  beam  onto  the  specified  target  sensor. 
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VIEWED  FROM  EARTH 


(»)  Acoujto-opiictlly  moduUtcd  tad  putted  tetuor  potkiocu 
Fig.  I  —  Schcmuic  drawing  of  LACE  teotor  inty 
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VIEWED  FROM  BARTH 

(b)  Mcehsaksliy-tnodutated  tabued  tenter  portions 
Fig.  1  {Continued)  —  Scherauie  drswtag  of  LAC £  sensor  imy 
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Fi*.2  —  Schematic  drawing  of  LATS  configuration  used  lor  the  k*eg ration  tests 
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Generation  of  the  required  radiation  was  accomplished  with  five  lasers:  25  Hz  excimer  and  10  Hz 
Nd:YAG  lasers  were  the  pulsed-radiation  sources;  25  kHz  acousto-optically  modulated  argon  ion, 
1440  Hz  mechanically  chopped  helium  neon  (HeNe),  and  1440  Hz  mechanically  chopped  deuterium 
fluoride  (DF)  lasers  were  the  continuous-wave  sources.  The  LACE  test  plan  [1]  describes  the  radiation 
specifications  and  Tables  1  and  2  summarize  them. 


Table  1  —  Laser  Characteristics 


Laser 

Wavelength 

_  (m)  m 

Beam  Diameter 
(mm) 

Maximum 

Power 

PRF 

(Hz) 

Pulse  Length 
(ns) 

Excimer 

0.308 

22 

90  mJ 

25 

19 

Ar+ 

0.5145 

3 

250  mW 

25k 

— 

Nd:YAG 

0,355 

8 

250  mJ 

10 

10 

Nd:YAG 

1.064 

8 

1  J 

10 

10 

DF 

3.6-4.2 

5 

1  W 

1440 

HeNe 

3.39 

5 

5  mW 

1440 

Table  2  -  LATS  Intel 

^ration  Requirements 

Laser 

Wavelength  (jsm) 

Power  Density  (W/cm2) 

Maximum 

Minimum 

Excimer 

0.308 

10* 

10? 

DF 

3.64.2 

10"4 

5  x  10-* 

Ar+ 

0.5145 

10~* 

O 

i 

© 

Nd:YAG 

0.355 

10* 

10? 

Nd:YAG 

1.064 

105 

10* 

The  goniometer  was  a  Klinger  BG  1 20  YZ  rotary  cradle  mounted  on  an  RT  120  XZ  rotary  stage  with 
a  Klingei  CC-1  programmable  indexer  and  a  step  resolution  of  0.01s.  Both  the  stage  and  cradle  had  an 
incremental  encoder  for  precise  closed-loop  positioning  and  speed  control.  These  two  rotary  stages  are 
combined  to  produce  independently  controlled  motion  in  azimuth  and  elevation.  The  goniometer  mirror 
was  in  a  two-axis  gimbie  mount  positioned  with  the  mirror  center  at  the  intersection  of  the  goniometric 
axes.  Two  independent  stepper-motor  channels  controlled  the  motion  of  the  two  orthogonal  axes  and 
allowed  for  manual  as  well  as  programmed  positioning.  The  stepper-motor  controller  is  able  to  establish 
a  home  reference  position,  make  absolute  movements  relative  to  that  position,  or  move  incrementally 
relative  to  a  given  position. 

For  LATS,  timing  control  was  provided  by  an  electronic  circuit  rather  than  the  specified  LACE 
trigger  laser.  Sync-out  pulses  from  the  pulsed  lasers  and  motor -driven  choppers  provided  reference 
signals  to  trigger  sensor  panel  data  collection.  The  timing  circuit  could  function  at  any  of  the  necessary 
detection  frequencies.  No  timing  circuit  was  necessary  for  the  acousto-optically  chopped  laser  because 
of  the  short  modulation  times;  a  data  pulse  was  always  present  in  the  detection  window. 
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The  LATS  computer  was  an  IBM  PC-AT  interfaced  with  two  IEEE-488  buses.  IBM  general  purpose 
interface  bus  (GPIB)  adapters  were  installed  in  the  PC  for  each  required  interface.  Through  the  GPIB, 
the  computer  communicated  with  the  goniometer  controller,  the  A/D  converter,  and  the  LACE  microvax 
sensor  array  subsystem  test  unit  (SASTU)  as  Fig.  3  schematically  shows.  Two  independent 
communication  buses  were  used  in  the  PC  to  avoid  conflicts  in  the  working  environments  and  to  reduce 
processing  time.  The  first  bus  served  as  a  data  link  between  the  LATS  PC  and  the  SASTU.  The  second 
bus  linked  the  PC  to  the  stepping-motor  controller  and  the  Hewlett-Packard  593 13A  A/D  converter.  The 
SASTU  served  as  system  controller-in-charge  (CIC)  on  the  first  bus.  The  IBM  PC-AT  LATS  computer 
was  controller  on  the  second  bus  and  operated  the  two  independent  axes  of  the  goniometer  and  the  A/D 
converter  as  three  independent  devices.  Table  3  summarizes  the  contents  of  the  configuration  files 
(IBCONF)  for  the  devices  used  with  the  GPIB  bus.  Reference  2  provides  a  description  of  GPIB  control. 
Appendix  A  lists  the  required  subroutine  sequence  for  goniometer  movement,  data  acquisition,  and  data 
transfer. 


Fig.  3  — Schematic  drawing  of  LATS  system.  The  Uacr  bum  is  shown  in  bold. 


Table  3  —  GPIB  Configuration  Parameters  for  LATS 


Device 

Address 

Time  Out 
(s) 

EOS*  Byte 

CIC 

EOIt 

Access 

GPIBO 

00H 

3 

01H 

YES 

YES 

mm 

GPIB1 

01H 

none 

OAH 

NO 

YES 

hi 

XAXIS 

06H 

3 

01H 

YES 

GPIBO 

YAXIS 

07H 

3 

01H 

YES 

GPIBO 

ADC 

08H 

10 

0IH 

YES 

GPIBO 

MVAX 

3 

OAH 

YES 

GPIB1 

•End  of  suing 

tSend  EOi  (cad  of  data  message)  with  EOS 
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Programmed  control  consisted  of  an  ASCII-coded  set  of  high-level  instruction  commands  written  to 
the  goniometer  controller  from  the  LATS  computer  over  the  interface  bus.  These  ASCII  strings 
controlled  goniometer  device  velocity,  acceleration,  displacement,  directionality,  scanning  interval,  and 
read/write  information  to  the  devices.  The  CC-1  stepper-motor  controller  manual  gives  a  full  description 
of  the  programmed  and  manual  operation  of  the  stepper-motor  controller. 

3.  POINTING 

Initially,  it  was  thought  that  accurate  pointing  could  be  achieved  by  evaluating  the  simple  planar 
trigonometric  relations  governing  a  point  source  incident  on  a  flat  target.  The  computation  of  the  tangent 
of  the  enclosed  angle,  the  ratio  of  the  vertical  and  horizontal  distances  of  the  sensor  from  the  origin,  and 
the  distance  from  the  steering  mirror  to  the  target  panel  yielded  azimuth  and  elevation  angles 
corresponding  to  the  pointing  vector  for  a  given  sensor  location.  A  scaled-down  laboratory  demonstration 
revealed  the  large  inaccuracies  of  this  simple  approach. 

The  task  of  pointing  a  beam  steered  by  an  azimuth-elevation  device  onto  a  planar  surface  at  a  close 
distance  is  to  map  a  spherical  surface  onto  an  intersecting  plane.  Since  the  sensors  on  the  test  panel  were 
arranged  in  a  flat  plane  indexed  by  Cartesian  coordinates  and  the  motion  of  the  incident  light  swept  out 
spherical  contours,  a  tangential  geometric  distortion  was  introduced,  and  it  was  necessary  for  an 
appropriate  mapping  to  be  derived. 

Figure  4  shows  the  laboratory  coordinate  system  used  to  characterize  the  mapping  of  a  spherical 
surface  onto  a  planar  surface.  The  origin  is  at  the  center  of  the  steering  mirror,  coincident  with  the 
goniometer  origin;  angles  are  measured  in  the  right-handed  sense  and  the  {-direction  toward  the  target 
is  negative. 

♦ 


+  2 

V 


fig.  4  —  Ltbonlof?  coontoute  u*e«3  for  LATS 
For  notation  purposes,  the  generalized  vector  v  is  defined  in  terms  of  the  unit  vector  as 


vj  ♦  vj  ♦  Vtt  =  (v„vrvs)r 

(1) 

v  rn  V0  -  V<P,,  Pt)r, 

<2) 

where  (v,,  vy  v.)r  is  the  notation  for  the  transpose  of  the  column  vector  for  malrix  manipulations.  The 
magnitude  of  the  vector  v  is  given  by 
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In  general,  as  Fig.  5  shows,  a  reflection  of  a  mirror  is  described  with  three  vectors,  s  is  the  vector  from 
the  steering  mirror  to  the  source,  T  is  the  vector  from  the  steering  mirror  to  the  target,  and,  IT  is  the 
mirror  normal  vector  that  is  the  coplanar  bisector  of  the  angle  formed  by  ~s  and  ~r  shown  in  Fig.  S. 


Unit  vector  toward  target 
Hr,  Target 

^ >*»  \ 

■  n - W — ► 

^Derived  unit  normal 

O 

i  "s  Source 
Unit  vector  toward  source 

Fig.  5  —  Schematic  repreaenutioa  of  the 
vectora  for  LATS  point  tag  algorithm 

Two  steps  are  required  to  determine  the  transformation  to  map  from  the  goniometric  system  to  a 
planar  surface.  The  first  step  is  to  find  where  the  mirror  normal  d  should  be  located  to  reflect  light  off 
the  goniometer  mirror  from  the  source  at  S  to  a  target  at  t ,  and  the  second  step  is  to  find  the  angles  (4, 
6)  that  will  position  die  goniometer  steering  mirror  normal  at  this  /l. 

The  accuracy  of  die  pointing  algorithms  depends  on  the  accuracy  with  which  the  experiment 
components  are  located  and  oriented  in  the  laboratory  system.  For  the  goniometer,  this  is  implicit  and 
fixed.  By  design  its  coordinates  define  the  laboratory  system,  and  other  LATS  components  must  by 
calibrated  with  respect  to  it. 

Ideally,  the  target  face  is  perfectly  flat,  is  perpendicular  to  the  laboratory  negative  i-axis,  and  has 
x *  and  y~axes  parallel  to  the  laboratory  axes  defined  by  the  goniometer.  The  target  face  flatness  is  least 
amenable  to  calibration.  In  the  LATS,  the  target  face  Is  the  sensor  panel  of  a  precisely  machined  satellite 
that  is  large  and  massive  but  not  under  unreasonable  strain.  Hie  width  of  the  beam  relative  to  the  size 
of  the  irradiated  sensor  area  made  target  flatness  a  negligible  correction. 

Perpendicularity  of  the  target  face  and  the  negative  laboratory  z-ax  is  is  achieved  to  adequate  accuracy 
through  mechanical  placement  of  the  two  systems.  This  also  ensures  collinearity  of  the  target  and 
goniometer  origins.  However,  it  is  also  a  simple  matter  to  redefine  the  target  face  coordinate  system  to 
lie  on  the  laboratory  z-axis.  Then  the  new  target  sensor  positions  are  computed  from  the  new  target 
origin  offset. 


Finally,  the  target  face  x-  and  y-axes  must  be  parallel  to  the  laboratory  axes  and  not  rotated  by  angle 
a  about  the  z  axis.  Two  approaches  are  used  to  correct  this  error:  position  the  target  and  the  goniometer 
so  that  the  axes  are  parallel,  or  measure  and  correct  for  a  relative  displacement  by  rotating  the  desired 
target  face  {xp  coordinates  through  -*a  into  the  laboratory  frame  coordinates  (xL,  yj  and  pointing 
the  beam  at  (jcL,  yj  instead  of  yj  for  each  sensor: 


XL 

cos  or  -  sin  a 

XT 

m 

yL 

sin  a  cos  a 

» 

(♦) 


Either  approach  requires  calibrating  the  source  steering  mirror  as  discussed  in  the  LATS  calibration 
section  below. 
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To  measure  the  rotational  misalignment,  direct  the  goniometer  to  a  convenient  (jt,  y)  position  in  the 
laboratory  frame.  (0, 0)  is  suggested  to  simplify  the  expressions.  Rearranging  the  expression  in  Eq.  (4) 
yields  the  following  expressions  for  a.  At  x  =  0, 


XT 

sin  a  =  — , 

y 

where  y  is  the  measured  target  frame  coordinate,  and  at  y  =  0, 


(5) 


(6) 


Similarly,  x  is  measured  in  the  target  frame  coordinate  system.  The  beam  is  pointed  toward  large 
excursions  along  the  x-  and  y-axes,  and  deviations  xT  and  yT  from  the  targets  are  measured  and  used  to 
compute  a.  If  a  nonzero  a  is  determined,  target  frame  coordinates  are  transformed  to  the  laboratory* 
parallel  target  frame  by  using  Eq.  (4). 


Determination  of  the  Mirror  Normal 

The  desired  normal  orientation  of  the  mirror  normal  vector,  in  terms  of  components  along  the 
laboratory  axes,  is  given  as  the  conlanar  bisector  of  f  and  S , 


where 

«-(»•*>’•  (ia/2)<i  ♦  v,  ♦  y,  ♦ 

B-teause  the  source  beam  is  reflected  one  or  more  times  before  impinging  on  the  steering  mirror  at  the 
goniometer  origin  the  light  source  7  is  located  in  the  laboratory  frame  at  the  center  of  the  TM  as  Fig. 
2  indicates.  The  vector  ~r  is  defined  similarly  as  the  sensor  largei  position  in  the  laboratory  frame,  with 
the  origin  of  the  vector  r  at  the  center  of  the  steering  mirror  and  the  target  at  the  center  of  the  sensor. 
By  measuring  the  source  and  target  vectors,  the  terms  s,  S,  r,  and  t  are  easily  computed  from  Eqs,  (2) 
and  (3),  and  (he  mirror  normal  is  then  determined  from  Eqs.  (7)  and  (8). 

This  approach  accounts  for  arbitrary  misalignments  in  almost  all  angles  and  positions  with  the 
exception  of  two:  the  steering  mirror  must  be  at  the  origin  of  the  goniometer  so  that  the  vertex  of  the 
reflection  is  at  the  center  of  rotation,  and  the  transfer  mirror  is  presumed  to  reflect  the  light  beam  directly 
onto  the  origin.  Arbitrary  misalignments  are  accommodated  as  long  as  the  angles  required  to  reflect  the 
User  beam  from  the  sensor  array  origin  to  the  source  can  be  measured  accurately. 

Determination  of  the  Angles  to  Position  the  Mirror  Normal 

Finally,  the  angles  ($,  6)  that  will  position  A  for  the  chosen  target  must  be  determined.  These  angles 
will  be  dependent  on  where  the  mirror  normal  is  when  ($,  6)  —  (0,  0).  This  is  defined  as  the  initial 
orientation  mirror  normal.  %  The  goniometer  angles  («^,  ^  are  the  angles  that  reflect  the  laser  beam 
from  the  origin  directly  back  to  the  source  s  on  the  transfer  mirror.  Again,  4g  and  (4o*  V  ***  measured 
in  the  laboratory  frame. 
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The  standard  matrix  presentation  of  a  rotation  of  a  unit  vector  0  through  6  about  x  is 

1  0  0  Kj  % 

j^(5)0  =  0  cos 5  -sin6  ^  =  t^cosfi  -  t^sinfi  , 

0  sinfi  cosfi  t>:  0ysin5  ♦  Pzcosd 


Similarly,  for  a  rotation  through  #  about  y, 

cos#  0  sin# 

J$(#)t>  «  0  10  % 

-sin#  0  cos#  $ 


tfxcos#  ♦  t^cos# 


-t^sin#  ♦  t^cos# 


Note  that  the  elevation  rotation  leaves  x  unchanged,  and  the  azimuthal  rotation  leaves  y  unchanged,  as 
these  are  the  corresponding  orthogonal  rotation  axes.  The  combined  rotation  is  given  by  the  product 
matrix. 


£(#,6)0  -  /?(#)^(6)P 


cos#  sin#stnS  sin#cos6  {A 


-sin#  cos#sin6  cos#  cos  <5  (  p 


/?■'(#,  6)  “  1$  r(#,  6)  for  this  orthogonal  matrix,  thus  simplifying  the  matrix  transformations  required 
to  compute  the  mirror  normal. 


rt(#,S) 

V 

% 

et 

t^ous#  *  ^ sin# sin 6  ♦  P.sin#cos5 
<Lcosd  -  P.sinfi 

J  * 

-^sin#  ♦  PJ,cos#sin6  ♦  ^cos#cos6 

To  determine  what  angles  are  required  to  direct  the  source  beam  onto  the  target,  the  rotation  that 

carries  the  initial  mirror  normal  Aq  -  (d^  d^,  A^f  into  the  desired  minor  normal  d  -  (/!,,  Ay,  A.f 
is  determined  from  Eq.  (12): 


d,  *  dicers#  ♦  A,^sin#sin6  ♦  d0tsin#cosd 
^  -  Accost  -  dfeSinS 


A.  °  -  d^stn#  ♦  d^c0s#sia6  ♦  dQ.cos#cos5.  (15) 

The  trigonometric  functions  in  #  and  6  given  in  Eqs.  (13)  and  (15)  complicate  rearrangement  into  closed 
expressions  for  the  angles  themselves;  some  trigonometric  substitutions  and  the  Pythagorean  relation  are 
required  for  solution;  the  details  are  described  in  Appendix  B, 
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8  =  COS'1 


*> 

—  cos"1 

% 

(1  - 1&)\ 

a  -  4)’_ 

(16) 


<j>  =  cos'1 


*z 

-  cos-1 

sin  6  +  ^  cos  8) 

.(1  " 

(1  -  A?)' 

(17) 


When  the  goniometer  is  positioned  at  (0,  8),  light  from  the  transfer  mirror,  positioned  at  (sx,  sy,  sz)  and 
incident  on  a  mirror  at  the  goniometer  origin  at  orientation  (%.,  n^)  when  (<£,  8)  =  (0,  0),  is 

reflected  to  the  target  position  (rx,  ry,  rz). 

The  mirror  normal  vector  required  to  direct  the  source  onto  the  target  is  easily  computed  from  the 
values  of  the  source  and  target  normal  vectors 


2 


i  =xty,z 


(18) 


and  normalizing  the  resultant  vector. 

LATS  Calibration 

It  is  important  to  precisely  determine  steering  mirror  position  and  to  accommodate  imperfect 
goniometer  mirror  mounting.  The  following  procedure  does  not  depend  on  the  presumption  of  target 
flatness  nor  on  the  determination  and  elimination  of  the  rotational  alignment  error  a  that  are  both  zero 
at  the  origin.  It  is  an  iterative  computation  dependent  on  the  two  sets  of  goniometer  angles  (</>,  8)  that 
reflect  the  source  beam  onto  the  target  origin  and  the  steering  mirror  in  the  laboratory  frame. 

The  steering  mirror  is  placed  in  the  -z  hemisphere,  as  close  angularly  to  the  target  origin  as  possible 
without  eclipsing  a  sensor.  Since  the  goniometer  mirror  is  to  reflect  the  beam  at  twice  the  normal  angle 
of  incidence,  there  is  an  angular  precision  bias  between  targets  on  the  source  mirror  side  of  the  target 
origin  and  targets  on  the  opposite  side  of  the  target  origin.  Angular  inaccuracy  and  imprecision  are 
doubled  by  the  reflection  and  are  proportional  to  the  source-goniometer/target-angle  magnitude. 
Therefore,  all  angles  are  minimized  by  keeping  the  steering  mirror  angularly  close  to  the  target  origin. 

the  steering  mirror  cannot  be  located  angularly  coincident  with  the  target  origin,  the  alternatives  are 
to  locate  it  in  the  os-plane  eliminating  y  bias,  or  in  the  yc-plane  eliminating  x  bias.  Because  of  the 
noncritical  accuracy  requirements  of  the  present  experiment,  the  steering  mirror  was  placed  approximately 
8°  from  the  target  origin. 

Two  simple  measurements  on  the  LATS  are  sufficient  to  permit  an  iterative  computation  of  the 
mirror  normal  that  reflects  a  source  beam  onto  the  target  origin  at  (0, 0,  -z).  These  two  measurements 
are  the  goniometer  angles  required  to  reflect  the  source  beam  onto  the  target  origin,  and  the  angles 
required  to  reflect  the  souice  beam  back  onto  itself.  Despite  .^accuracies  in  initial  estimates  of  both  the 
source  position  and  the  at-rest  goniometer  mirror  normal  /Iq,  the  prescribed  measurements  contain 
sufficient  information  about  the  angle  to  converge  numerically  to  a  self-consistent  and  correct  value  of 
both  position  vectors  simultaneously. 
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It  is  convenient  but  not  important  that  (<j>,  S)  be  nearly  zero.  The  goniometer  used  in  the  LATS 
supported  a  logical  definition  of  (<t>,  6)  =  (0,  0)  that  was  used.  This  constitutes  the  initial  orientation  of 
the  goniometer  mirror  normal.  The  estimated  position  of  the  source  (x,  y,  z)  is  measured  in  the 
laboratory  frame.  The  vector  is  normalized  and  the  result  designated  S0  where  the  subscript  denotes  an 
initial  approximation  to  the  true  source  normal .  It  is  precisely  because  of  the  difficulty  of  measuring  both 
of  these  critical  parameters  accurately  that  the  method  described  is  used. 

The  first  iterative  correction  to  the  initial  mirror  normal  vector  pointing  at  the  sensor  array  origin, 
r  ss  (0,  0,  -1),  is  given  as 

=  (Sqx  -  0,  sQy  -  0,  -  1)  =  (2  -  2$$^  (sq,,  sQy,  Sqz  -  1)  (^) 

from  Eq.  (7),  with  the  rotation  given  by  Eq.  (12)  that  reflects  the  source  beam  onto  the  origin.  This 
new  Aq  is  used  in  Eqs.  (13)  and  (1.5)  to  determine  a  new  £0  from  the  measured  angles  (^,  Sq)  required 
to  reflect  the  initial  mirror  normal  onto  the  source.  These  values,  in  turn,  are  used  to  iteratively 
determine  a  new  initial  mirror  normal  in  Eq.  (19)  until  the  solution  converges. 

Because  of  the  angle-doubling  effect  of  beam  reflection,  an  error  in  the  mirror  normal  causes  twice 
the  error  in  the  reflected  beam.  Conversely,  an  error  in  the  reflected  beam’s  measurement  corresponds 
to  half  the  error  in  the  mirror  normal.  Beginning  the  iterative  procedure  with  the  right-hand-side  of  Eq. 
(19)  supports  convergence  bee-  ,  e  i  ors  in  the  position  of  $  are  halved  before  correcting  with  Eqs.  (13) 
and  (15)  on  the  next  iteration.  Convergence  is  guaranteed  if  the  experiment  is  designed  with  the 
goniometer  at-rest  mirror  normal  /Iq  .  not  the  source  position  1.  which  is  angularly  closer  to  the  laboratory 
target  origin.  Then  the  rotation  required  to  reflect  the  mirror  normal  onto  the  target  origin  is  smaller  than 
that  required  to  reflect  the  mirror  normal  onto  the  source.  For  the  nominal  position  /Iq  this  error 
is  halved,  which  roughly  cancels  the  angle-doubling  error  introduced  in  Eqs.  (13)  and  (15).  Application 
of  Eq.  (19)  halves  error,  while  application  of  Eqs.  (13)  and  (15)  preserves  them,  and  the  procedure 
converges. 

Implementation  of  this  goniometric  mapping  in  the  laboratory  demonstrated  a  large  improvement  over 
the  initial  approach.  However,  the  accuracy  of  the  mapping  is  a  function  of  the  accuracy  of  the  input 
parameters,  including  presumptions  of  perpendicularity  of  the  target  plane  and  the  laboratory  z-axis.  For 
the  supplied,  not  measured,  sensor  locations  used  here,  each  position  was  indexed  by  its  xy  distance  from 
the  satellite  origin  and  specified  to  accuracies  of  1.0  cm  for  the  infrared  sensors  and  0.5  cm  for  the 
ultraviolet  and  visible  sensors;  errors  larger  than  those  specified  are  thought  to  occur  for  several  sensor 
locations  in  the  model  spacecraft  used  for  the  integration  tests.  These  unexpectedly  large  inaccuracies 
in  sensor  placement  from  the  indexed  Cartesian  coordinates  required  the  mapping  to  be  augmented  with 
manual  position  fine  tuning  to  establish  the  final  goniometric  pointing  vector. 


4.  POWER  COMPUTATION 

Incident  power  computations  were  made  by  the  LATS  system  and  sent  to  the  SASTU  for  comparison 
with  the  powers  collected  at  the  sensors.  The  ratio  of  the  beam  power  reflected  by  the  beamsplitter  to 
that  transmitted  through  the  optical  path  to  the  target  was  known  with  accuracies  of  10%  for  each 
wavelength  used.  The  computation  of  incident  power  required  knowledge  of  laser  and  sensor-dependent 
parameters.  These  parameters  included  any  neutral  density  filtering  in  the  laser  beam  (ND),  the  pulse 
repetition  frequency  (PRF)  of  the  laser,  the  Dulse  length  (PL)  of  the  laser,  and  the  ratio  of  the  areas  of 
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the  sensor  and  the  laser  beam.  Equations.  (20)  and  (21)  give  the  incident  power  for  the  chopped  and 
pulsed  sensors,  respectively: 


Wi  ■  A.  •  10'“  • 


V*  *  A,  '  10*  •  PRF-'  '  PL 


Ar^unscr 

Area^  ’ 

(20) 

-I  .  sauor 

(21) 

Pm  is  the  calibrated  power  reading  from  the  analog  to  digital  converter.  To  optimize  the  placement  of 
the  laser  beam  at  each  sensor  following  goniometric  and  manual  pointing,  the  beam  was  scanned  over 
a  0.25°  x  0.25°  grid  in  5  x  5  steps  centered  on  the  sensor.  Average  power  levels  are  measured,  and 
the  position  yielding  the  highest  power  is  retained  for  the  final  vector. 

5.  LATS  ALGORITHM 

Three  computer  programs  were  used  for  the  LATS  tests.  An  initialization  routine  created  and 
formatted  the  five  sensor-dependent  position  information  bases  for  the  335  satellite  sensors  into  five 
sensor-dependent  location  files  (pulsed  high  power,  pulsed  low  power,  modulated,  infrared  high  power, 
and  infrared  low  power).  The  infrared  and  ultraviolet  sensors  are  distinct  detector  pairs  that  cover  well- 
defined  power  ranges.  A  setup  routine  calculated  the  goniometer  pointing  vectors,  allowed  manual  vector 
correction,  and  wrote  the  final  vector  database.  Finally,  a  pointing  routine  interfaced  with  the  SASTU 
passed  sensor  identification  numbers,  pointed  the  laser  beam  at  the  specified  sensor,  and  passed  power 
data  back  to  the  SASTU.  All  programs  were  written  in  RM  Professional  Fortran  (3],  and  the  setup  and 
pointing  routines  are  included  in  Appendix  C. 

Initialization  Routine 

The  following  describes  the  procedures  in  the  initialization  routine. 

A.  Prompt  user  for  input  laser  identification  number  to  distinguish  among  the  five  sensor-dependent 
location  files. 

B.  Input  Cartesian  coordinates  of  the  sensor  with  respect  to  the  satellite  origin. 

C.  Write  the  indexed  coordinates  to  the  sensor-dependent  location  file. 

Setup  Routine 

The  following  describes  the  procedures  in  the  goniometer  setup  routine. 

A.  Prompt  user  for  the  sensor  identification  number  and  laser  type. 

B.  Look  up  indexed  position  information  from  sensor  location  file. 

C.  Compute  and  calibrate  pointing  angles. 

1 .  Compute  rest  mirror  normal  /Iq  from  steering  mirror  and  target  vectors. 

2.  Compute  desired  normal  orientation  ft  from  Eq.  (7)  knowing  the  target  vector  and  the 
source  vector.  The  initial  mirror  normal  is  computed  by  using  the  solution  to  recoraput^f 
that  is  then  solved  again  for  the  mirror  normal  until  the  solutions  converge. 

3.  Compute  elevation  and  azimuth  angles  for  sensor  location  from  Lqs  (16)  •md  (!7). 
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D.  Direr*  the  laser  beam  onto  the  specified  sensor. 

1.  Covert  angular  information  into  step  interval  and  form  ASCII  string  commands.* 

2.  Writ^  the  command  string  to  the  goniometer. 

E.  Prompt  user  to  manually  correct  laser  beam  placement  on  the  sensor. 

F.  Scan  across  sensor  face  to  locate  laser  beam  placement  for  maximum  sensor  irradiance  (optional 
routine). 

G.  Record  position  of  maximum  irradiance  on  position  vector  database. 

Pointing  Routine 

The  following  describes  the  procedures  followed  in  the  goniometer  pointing  routine. 

A.  Read  sensor  identification  from  SASTU. 

B.  Obtain  sensor  azimuth  and  elevation  angles  from  vector  information  database. 

C.  Direct  laser  beam  r  n  the  ;pecified  sensor. 

D.  Recokd  incident  laser  beam  power. 

E.  Loop  to  A. 

Upon  completion  of  the  ..utomated  individual  sensor  tests,  a  flood  test  was  performed  at  each 
wavelength.  3y  reflecting  the  laser  beam  from  the  transfer  mirror  to  the  convex  mirror/concave  lens 
assembly  illustrated  in  Fig.  2,  the  beam  wus  expanded  to  irradiate  all  sensors  on  the  satellite 
simultaneously  for  the  multisensor  flood  tests.  These  rests,  used  to  determine  whether  there  is 
interference  between  sensors,  do  not  require  precision  pointing  or  computer  control. 

6.  VACUUM  TESTS 

Figure  6  shows  several  modifications  that  were  made  on  the  lATS  for  the  vacuum  chamber  test?. 
A  doubled  Nd:  YAG  laser,  an  argon-ion  laser ,  and  an  infrared  quartz-halogen  lamp  were  the  three  sources 
used  for  the  vacuum  tests.  Because  of  mechanical  limitations,  the  infrared  lamn  was  used  only  for  flood 
testing  of  the  sensor  array,  and  the  two  visible  lasers  were  used  for  both  multiscnsor  flood  testing  and 
individual  sensor  irradiation.  Interchangeable  lenses  placed  he  entrance  port  to  the  vacuum  chamber 
provided  the  collimated  or  diverging  laser  beam  for  the  two  types  of  tests,  individual  and  flood. 

Long  periods  of  time  were  required  fr  attain  the  vacuum  in  the  large  satellite  chamber.  To  avoid 
having  to  open  the  chamber  between  thermal  cycles  ot  between  individual/flood  tests  modes  to  reestablish 
the  vacuum,  a  remotely  controlled  method  of  switching  between  multi-sensor  flood  testing  and  Individual 
sensor  testing  was  devised.  A  convex  mirror  mounted  to  the  backside  of  the  goniometer  gimble  mount 
allowed  one  hemisphere  of  the  goniometer  to  reflect  collimafod  laser  light  and  the  other  hemisphere  to 
reflect  rapidly  diverging  light.  To  switch  between  the  two  types  of  tests,  the  goniometer  carousel  had 


•To  write  t  tiring  to  the  goniometer  controller  the  length  of  the  string  mud  be  specified.  b  it  ncctUMy  for  (hit  program,  n 
written  in  RM  Professions]  Fortran,  *o  branch  into  individual  subroutine*  to  ■pnommndttr  the  utringtize  of  the  g^ityrK*^ 
oomnand. 
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only  to  be  rotated  by  180°.  The  infrared  lamp  was  in  a  pressurized  canister  and  positioned  in  the  vacuum 
chamber  to  provide  complete  irradiation  of  the  satellite  face.  Except  for  the  modifications  just  described, 
validation  tests  for  the  scanned  individual  sensors  and  multisensor  tests  equivalent  to  those  outside  the 
chamber  were  performed  on  the  spacecraft  in  the  vacuum  chandler. 

7.  CONCLUSIONS 

The  LATS  system  provided  precise,  repeatable,  and  automated  goniometrlc  vectoring  for  preflight 
system  integration  and  operational  evaluation  of  the  LACE  spacecraft.  The  pointing  vector  algorithm 
developed  for  the  gonio'  *eter  provided  placement  of  a  laser  beam  at  any  sensor  location  on  the  satellite, 
based  only  on  the  indexed  Cartesian  coordinates  of  the  sensor  and  the  distance  of  the  goniometer  steering 
mirror  and  transfer  mirror  from  the  satellite  sensor  plane  origin.  Large  error;  in  sensor  location  from 
the  indexed  coordinates  were  corrected  with  manual  vector  positioning.  Cnee  pointing  vectors  wore 
established  for  all  satellite  sensor  locatio  ns,  automated  laser  testing  was  performed. 
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Appendix  A 

GPIB  ROUTINES  FOR  LATS 


GPIB  Routine  for  LATS  Goniometer  Motion 

C  move  x-axis  n  steps  (n  =  two  digit  integer,  e.g.  10) 
integer*2  xaxis 
character's  X$ 

C  assign  unit  descriptor  for  x-axis 
xaxis  =  ibfind  (’xaxis’) 

C  clear  the  device 
call  ibclr  (xaxis) 

C  file  'RATE"  contains  goniometer  velocity  and  acceleration  values 
C  write  file  "RATE"  to  the  x-axis 
call  ibwrtf  (xaxis,  rate) 

C  clear  the  device 
call  ibclr  (xaxis) 

C  construct  string  containing  N  number  of  steps 
X$  ■  ’N’//’  7/’10’//char(13) 

C  write  X$  to  x-axis 
call  ibwrt  (xaxis,  X$,  S) 

GPIB  Routine  to  Read  the  Analog  to  Digital  Converter 

C  read  converted  measurement  from  HP  converter 
integer*2  adc 
character *2  buffers 

C  assign  unit  descriptor  for  analog  to  digital  converter 
adc  =  ibfind  Cade’) 

C  command  for  converter  to  write  data  is  ’H8AJ' 
call  ibwrt  (adc,  ’H8AJ\  4) 

C  data  is  written  as  binary-coded  decimal  in  ASCII 
call  ibrd(adc,  buffers,  2) 

GPIB  Routine  to  Write  and  Read  to/from  the  SACTU 

C  assign  unit  descriptor  for  adapter  linked  to  SASTU 
integer  *2  gp 

character*!!)  bufferb,  bufferc 
gp  a  ibfind  CGP1B2’) 

C  write  data  to  GPIB  adapter 
call  ibwrt  (gp,  bufferb,  10) 

C  read  data  from  the  GPIB  adapter 
call  ihrd  (gp,  bufferc,  10) 
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DETAILED  SOLUTION  FOR  EQUATIONS  (16)  AND  (17) 


To  solve  for  0  and  5,  three  equations  are  used  in  two  unknowns.  These  trigonometric  representations 
must  be  arranged  into  closed  expressions  for  the  angles.  By  using  trigonometry  and  the  generalized 
Pythagorean  relation 


these  expressions  can  be  found. 
From  Eq.  (14), 


/ly  -  A*  cos  5  -  /Iq.  sin  6. 


(B2) 


By  using  the  trigonometric  construction  shown  in  Fig.  BI  and  making  the  following  substitutions, 
a  »  4^,  b  *  /Iqu  yields 

(B3) 


by  trigonometric  identity; 


*  c  (cos  a  cos  6  -  sin  a  sin  fi) 
-  c  •  cos  (a  ♦  6) 

■  (a1  ♦  62)1  cos  (or  *  6) 

Ay  •  toj  ♦  4>’cos(a  ♦  6). 


CM) 

(B5) 

(B6) 

(B7) 


c  ■  (a2 ♦  b  ) 
cos  a  «  a/c 
sina  *  h/c 


Pit-  Bl  —  Trigonometric  coaUruction  needed  to  timplify  the 
expreukxu  uied  to  eotve  for  Eq*.  (16)  tad  (17) 
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By  using  the  identity  in  Eq.  (Bl)  and  rearranging: 


cos  (a  +  4) 


(1  -/&)’ 


8  =  cos"1 


a-4)’ 


-  a 


8  -  cos"1 


a-4)’ 


-  cos 


-1 


V 


a  -4)’ 


m 

(B9) 

(BIO) 


This  is  a  closed-form  expression  for  6. 

From  Eq.  (15), 

At  =  m$  +  AQyCO$$s\n8+Af3zcas4>GQs6  (Bl  1) 

-  cos  ^(/^  sin  8  *  Aq.  cosd)  -  sin  M^),  (B12) 

which  is  of  the  same  form  as  Eq.  (B2).  Letting  a  “  sin  6  ♦  cos  6)  and  Z?  »  rearranging 
the  expression  yields 


A.  -  [(/loy  sin  6  ♦  Aq.  cos  8?  ♦  4]*  605  (4  ♦  ft- 
Simplification  is  possible  by  expanding  the  square-root  term  as  follows: 

(/toy  sin  6  ♦  Aq.  cos  4)2  ♦  4  *  4  4  *  4  0052  *  *  sin  4  «m  8  4 

-  4(1  -  cos2  8)  ♦  4(1  -  sin2  4)  ♦  4  +  2/V*ct  sin  4  cos  5 

-  (4  ♦  4  ♦  4)  ♦  sin  4  cos  4  -  4  6051  4  -  4  *in2  4 

-  1  -  (4  «>s2  4  -  24^  sin  4  cos  4  ♦  4  sin2  4) 

“  1  -  (4  cos  4  -  4^  sin  4)1 

-1-4 


0*13) 

0*14) 

0*15) 

0*16) 

(BIT) 

0*18) 

0*19) 
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from  Eq.  (B2).  From  Eq.  (B13) 


*  (1  -  /lyV  cos(*  ♦  0) 


(B20) 


$  *  cos 


-l 


•  • 

z 

-  cos"1 

0  -  «,v 

(/^  sin  5  ♦  ^  cos  5) 

(i  -  «?>’ 


(B21) 


Eq.  (B21)  is  a  closed  form  expression  for  4>. 


n  n 


APPENDIX  C 


LATS  PROGRAMS 


Cl.  Setup  Routine 
C  19  NOV  1937 

C  PROGRAM  TARGETS ET.FQR  SEPT1 , 1987 

C  MODIFIED  FROM  PROGRAM  LATS. FOR  AUGUST,  1987 

C  MUST  BE  LINKED  WITH  FILES  PFIB,  M1SET,  M2SET,  SCAN 
C  PROGRAM  SERVES  AS  AN  INTERACTIVE  SET  UP  PROGRAM 
C  REQUESTING  SENSOR  #,  MOVING  GONIOMETER  TO  SENSOR  LOCUS 
C  PROMPTING  USER  TO  MANUALLY  POSITION  BEAM,  SCANNING  BEAM 
C  ACROSS  SENSOR  TO  MAXIMIZE  1RRADIANCB,  AND  RECORDING 
C  THIS  POSITION  INFORMATION  ON  DATA  FILES 

C  THE  FOLLOWING  IS  GPIB  PROFESSIONAL  FORTRAN  HEADER 
COMMON  /IBGLOB/  IBSTA 
integers  ibclr,  TIMO,  END 
integer*!  ibftnd,  IBSTA 
integer*!  ibrd,  ibrdf,  ibtmo 
integer*2  ibwrt,  ibwrtf 

INTEGER *2  TNONE,  T10us,T30us,T100us,T300us 
INTEGER*!  Tlms,T3ms,T10ms,T3ftns,T100ms 
INTEGER*!  T300ms,Tls,T3s,T10s,T30s 
INTEGER *2  T100s,T300s,T1000s,CMPL 
C  VARIABLE  DECLARATIONS  FOR  PROGRAM 
INTEGER  XAXIS,  YAXIS,  GP,  DN 
INTEGER*2  BD,  GPIBO,  IOUT,  F,  ENTRY 
INTEGERS  PPR,  SPR,  IXX,  IYY,  IX,  IY,  PI,  DI 
REAL  R(3),  SH(3),  MH(3) 

REAL  PH,  DL,  XP,  YP 
CHARACTER*!  SE,  S$t  TS,  US,  V$,  GS 
LOGICAL  BRANCH,  NEW 
c  EOS  mode:  bit  values 

data  BIN/ZM0007,XEOSyZW/,R£OS/Z’4007 
c  Timeout  values 

data  TNONE/0/,TlOus/l/,T3Ous/2/,Tl(X)us/3/,T3OOus/4/ 
dau  Tlras/5/,T3ms/6/,T10ms/7/,T30ms/8/,T100ms/9/ 
data  T300ras/10/,TIs/l  l/,T3s/12/,TI0s/13/,T30s/14/ 
data  T100s/15/,T300s/16/,T1000s/17/ 
data  T100s/lS/,LF/10/ 

GIVE  ONLINE  INSTRUCTIONS  TO  USER 

WRITE  (*,*)'  TO  RUN,  THE  FOLLOWING  FILES  MUST  BE  SUPPLIED:* 
WRITE  (*,*)'  "HARDWARE.DAT*  CONTAINS  SYSTEM  MEASUREMENTS:* 
WRITE  (*,*)*  FACTOR  (1.0),’ 
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WRITE  (*,*)’  HORIZ  DISTANCE  TO  BEAM  ORIGIN  FROM  PANEL  CENTER,’ 
WRITE  (*,*)’  VERT  DISTANCE  TO  BEAM  ORIGIN  FROM  PANEL  CENTER,* 
WRITE  (*,*)’  HORIZ  DISTANCE,  STEERING  TO  TRANSFER  MIRROR,’ 

WRITE  (*,*)’  VERT  DISTANCE,  STEERING  TO  TRANSFER  MIRROR,* 

WRITE  (*,*)’  Z  DISTANCE,  STEERING  TO  TRANSFER  MIRROR,’ 

WRITE  (*,*)’  ALL  DISTANCES  IN  CM  UNITS’ 

WRITE  (*,*)*  "DISPLAY"  CONTAINS  THE  COMMAND  "V  0[CR]*’ 

WRITE  (*,*)’  "RATE.DAT*  CONTAINS  A,R,S,F  GONIOMETER  PARAMETERS:’ 
WRITE  (*,*)’  "A[CR]R  2S0[CR]S  2[CR]F  20[CR]"’ 

OPEN  (20,  FILE=  ’DECQDE.DAT’,  STATUS = ’OLD’) 

CALL  SETUP  (F,R,SH,MH,XP,YP) 

XAXIS=  IBFIND  (’X AXIS’) 

YAXIS =  IBFIND  (’YAXIS*) 

GP=  IBFIND  (’GPIB2’) 

C  ASK  USER  TO  SELECT  DATA  ENTRY  MODE 
WRITE  (*,  *)  ’ 

WRITE  (*,  *)  ’  ENTER  (I)  FOR  KEYBOARD,  (2)  FOR  IEEE  DATA  BUS’ 

READ  (*,  *)  ENTRY 

C  PERFORM  DEVICE  CONNECTION  CHECK 

WRITE  (*,•)’  THE  FOLLOWING  TWO  RESPONSES  SHOULD  BE  "48"* 

WRITE  (*,*)’  IF  THEY  ARE  OTHERWISE,  POWER  DOWN  SYSTEM’ 

CALL  IBRSP  (XAXIS,  SPR) 

WRITE  (•,  •)  SPR 
IF  (SPR.NE.48)  GO  TO  1000 
CALL  IBRSP  (YAXIS,  SPR) 

WRITE  (*,  •)  SPR 
IF  (SPR.NE.48)  GO  TO  1000 

C  LOOP  HERE 

1  CALL  LOC  (F,IX,JY,DN,R,SH,MH,BRANCH,XP,YP,SE, ENTRY, GP) 

IF  (IX  .LT.  0)  GO  TO  100 
GOTO  UO 

100  S$=’-7/CHAR(!3) 

TS*’+*  //CHAR(I3) 

GO  TO  150 

110  S$!=*+7/CHAR(13) 

T$»’-7/CHAR(13) 

150  IF  (IY  .LT.  0)  GO  TO  200 

GO  TO  210 

200  US=’-7/CHAR(13) 

VS  *  ’  -f  7/CH  AR(1 3) 

GO  TO  250 

210  U$as'-f ’//CHAR(13) 

V$*’-7/CHAR(13) 

250  G$=’G7/CHAR(13) 

IF  (BRANCH)  GOTO  6 

CALL  MOVE2  (IX,  IY,  XAXIS,  YAXIS,  SS,  TS,  US,  VS,  GS) 

GOTO  7 

6  CALL  MOVE1  (IX,  IY,  XAXIS,  YAXIS,  SS,  T$,  US,  VS,  GS) 

7  CALL  POS  (SE,DN,XAXIS,YAXIS,GP,IX,IY .ENTRY) 

C  GO  TO  HOME 

IF  (BRANCH)  GOTO  8 
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CALL  MOVE2  (IX,  IY,  XAXIS,  YAXIS,  T$,  S$,  V$,  US,  G$) 

GOTO  9 

8  CALL  MOVE1  (IX,  IY,  XAXIS,  YAXIS,  T$,  SS,  VS,  US,  GS) 

9  IF  (ENTRY.EQ.3)  GOTO  10 
GOTO  1 

1000  WRITE  (*,  1001) 

1001  FORMAT  C  DEVICE  CHECK  FAILURE  ’) 

3  FORMAT  (3A) 

5  FORMAT  (LI) 

10  STOP 
END 

^  ************************************************************** 

Q  ************************************************************** 

q  ************************************************************** 

SUBROUTINE  LOC  (F,IX,IY,DN,R,SH,MH,BRANCH»XP,YP,SE,ENTRY,GP) 

C  ROUTINE  TAKES  MEASURED  GONIOMETER  PARAMETER  AND  CALCULATES 
C  PHI,  DELTA  FOR  GIVEN  SENSOR 

REAL  CD,  SD,  NN,  RR,  D 
REAL  R(3),  N(3),  RH(3),  SH(3),  NH(3),  MH(3) 

REAL  X,  XI,  X2,  A,  B,  DL,  PH 
REAL  XW,  YW,  P,  XP,  YP 
INTEGERS  DN,  I,  F,  IX,  IY,  SE 
INTEGERS  ENTRY,  DUM,  IAX,  1AY 
INTEGER *2  IBRD,  IBWRT,  IBFIND 
INTEGER  GP 
CHARACTER*:  UV,  IR 
CHARACTER*4  MESSAGE 
CHARACTER*6  FNAME.  POWER 
CHARACTER*  1 1  LIST 
LOGICAL  BRANCH 
IF  (ENTRY.EQ.l)  THEN 

C  ASK  WHICH  TYPE  OF  SENSOR  IS  SOUGHT 
SO  WRITE  (*,  5) 

READ  (*,  6)  SE 
WRITE  (*,  107) 

READ  (*,  *)  DN 
END  IF 

IF  (ENTRY.EQ.2)  THEN 
CALL  IBRD  (GP,  POWER,  6) 

WRITE  (*,  •)  •  DATA  READ  FROM  BUS  V 
CALL  DECODE  (POWER,  SE,  IAX,  IAY,  DN) 

WRITE  (*,  •)  SE,  DN 
MESSAGE  =>  ’INIT 

CALL  ENCODE  (IAX,  IAY,  SE,  MESSAGE,  UST) 

WRITE  (*,  •)  UST 
CALL  IBWRT  (GP,  LIST,  11) 

END  IF 
CLOSE  (7) 

IF  (SE.LE J.AND.SE.GT.O)  THEN 
FNAME=  ’UV.DAT* 

ELSE  IF  (SE.EQ.4.0R.SE.EQ.5)  THEN 
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FNAME=  TR.DAT 
ELSE 
GO  TO  50 
END  IF 

OPEN  (7,  FILE=  FNAME,  STATUS  =  ’OLD’) 

READ  (7,*)  I,  YW,  XW 
IF  (DN.EQ.I)  GO  TO  53 
GO  TO  52 
CLOSE  (7) 

R(l)=  XW  /  F 
R(l)=  R(l)  -  XP 
R(2)=  YW  /  F 
R(2)=  R(2)  -  YP 

RR=  SQRT  (R(l)**2  +  R(2)**2  +  R(3)**2) 

DO  551=1,3 
RH(I)=  R(I)  /  RR 
CONTINUE 

COMPUTE  REQUIRED  GONIOMETER  MIRROR  NORMAL,  EQS  1-3 
DO  60 1=1,3 
N(I)=  (RH(I)  +  SH(I))  12 
CONTINUE 

NN=  SQRT(N(1)**2  +  N(2)**2  +  N(3)**2) 

DO  65  1=1,3 
NH(I)=  N(I)  /  NN 
CONTINUE 

COMPUTE  DELTA,  EQN  19 
D=  SQRT  (1  -  MH(1)**2) 

X«  NH(2)  /  D 
XI  -  ACOS(X) 

X  =  MH(2)  /  D 
X2=  ACOS(X) 

DL=  X2  -XI 
SD=  SIN  (DL) 

CD=  COS  (DL) 

DL=  DL  •  ISO  /  3.14159 

C  COMPUTE  PHI,  EQN  24B 
A=  MH(1) 

B=  MH(2)*SD  +  MH(3)*CD 

XI  =  B*NH(1)  -  A*NH(3) 

X2=  A*NH(1)  +  B*NH(3) 

P=  (XI  /  X2) 

PH=  ATAN  (XI  fX2) 

C  CHANGE  SIGN  OF  PHI,  GONIOMETER  USES  POS  VALUES  FOR  CW 
PH=  -PH  •  180/3.14159 

C  CONVERT  TO  GONIOMETER  COUNTS  WITH  ROUNDING 
WRITE  (•,  111)  PH,  DL 
DL=  DL  •  100.0 
PH=  PH  •  100.0 
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IY=  ANINT(DL) 

IX*  ANINT(PH) 

IF  (IX.GT.-I00. AND.IX.LT.  100)  BRANCH  =  .TRUE. 

IF  (IX.GE.100.0R.IX.LE.-100)  BRANCH  =  .FALSE. 

WRITE  (*,  110)  IX,  IY 
FORMAT  (A) 

FORMAT  (’  IS  SENSOR  AR(1),UVH(2),UVL(3),IRH(4),IRL(5)?  ’/) 

FORMAT  Ul) 

107  FORMAT  (’  DETECTOR  NUMBER:  ’) 

1 10  FORMAT  C  NUMBER  OF  X  STEPS  ’,15  ,  ’  NUMBER  OF  Y  STEPS  \  15) 

1 1 1  FORMAT  C  AZIMUTH  PEG)’,  F6.2  ,  ’  ELEVATION  PEG)  \  F6.2) 
RETURN 

END 

Q  ft************************************************************ 

f  ********«********************«***«*****tt***3****************** 

SUBROUTINE  SETUP  (F,R,SH,MH,XP,YP) 

C  ROUTINE  GETS  SCALING  FACTOR,  GONIOMETER  POSTIONS  REQUIRED 
C  FOR  COMPUTATIONS 
INTEGER*2  F 
INTEGERS  RL 

REAL  DTR,  PH,  DL,  CD,  SD,  CP 
REAL  SP,  SS,  MM,  RR,  NN,  XP,  YP 
REAL  S(3),  M(3),  SH(3),  R{3) 

REAL  RH(3),  N(3).  NH(3),  MH(3) 

RL*  30 

OPEN  (2  .FILE = ’OAR . D AT’ , ACCESS  * 'DIRECT’.RECL « RL, 

/  FORM  » ‘FORMATTED*,  STATUS »  ‘OLD’) 

OPEN  (3 ,  FILE  *  ‘U  VH .  DAT’ .ACCESS  *  *DI RECT'  ,RECL  *  RL, 

/  FORM = ’FORMATTED*.  STATUS  =  *OLD‘) 

OPEN  (4,FILE=*UVL. DAT*, ACCESS 35 ‘DIRECT* ,RECL=RL, 

I  FORM « ‘FORMATTED’,  STATUS  «  ‘OLD*) 

OPEN  (8,  FILE *  *IRH .  DAT*  .ACCESS  « ‘DIRECT*  ,RECL = RL, 

/  FORM  « ‘FORMATTED*.  STATUS  «  *OLD*) 

OPEN  (9, FILE =TRL.DAT',  ACCESS*  ’DIRECT, RECL*  RL, 

/  FORM = ‘FORMATTED* ,  STATUS  =  *OLD*) 

OPEN  (10,  FlLEs  ’HARDWARE.DAT’,  STATUS  *  *OLD*) 

READ  (10,  *)  F,  XP,  YP,  S(l),  S(2),  S(3),  PH,  DL 
C  ITERATIVE  MIRROR  NORMAL  SOLUTION 

C  RUNS  500  ITERATIONS 

DTR*  3.14159/180 

C  INPUT  STEERING  MIRROR  COORDINATES 

C  INPUT  GONIOMETER  ANGLES  AT  R = (0,  0,  - 1 ) 

C  CONVERT  TO  RADIANS.  COMPUTE  TRIG  FUNCTIONS 

DL*  dl*dt; 

PH*  PH  •  DTR 
CD*  COS  PL) 

SD*  SIN  pL) 

CP*  COS  (PH) 

SP*  SIN  (PH) 
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C  NORMALIZE  STEERING  MIRROR  VICTOR 
SS=  SQRT(S(1)**2  +  S(2)**2  +  S(3)**2) 

DO  iO  1=1,3 

s(0=  SO)  /  ss 

10  CONTINUE 

DO  25  J= 1,500 

C  CALCULATE  GONIOMETER  MIRROR  NORMAL 

MM=  SQRT  (2*  (1-S(3))) 

M(l)=  (S(1)*CP  -  (S(3)-1)*SP)/MM 

M(2)=  (S(1)*SP*$D  +  S(2)*CD  +  (S(3)-1)*CP*SD)/MM 

M(3)=  (S(1)*SP*SD  -  S(2)*SD  +  (S(3)-1)*CP*CD)/MM 

C  COMPUTE  NEW  S  FROM  CALCULATED  GONIOMETER  MIRROR  NORMAL 
C  AND  MEASURED  PHI,  DELTA 

S(l)=  M(l)  *  CP  +  M(2)  *  SP  *  SD  +  M(3)  *SP*  CD 
S(2)=  M(2'»  *  CD  -  M(3)  *  SD 

S(3)=-M(l)  *  SP  +  M(2)  *  CP  *  SD  +  M(3)  *  CP  *  CD 
25  CONTINUE 
C 

C  COMPUTE  REST  MIRROR  NORMAL  NO  FROM  STEERING  MIRROR 
C  AND  TARGET  VECTORS 

SS=SQRT  (S(l)**2  +  S(2)**2  +  S(3)**2) 

DO  301=1,3 
SH(I)=  s©  /  ss 
30  CONTINUE 

C  INPUT  AND  NORMALIZE  TARGET  VECTOR 
R(l)=  0.0 
Rt?.)=  0.0 
k(3)  ~  -1.0 
RR  -  1.0 

C  READ  (*,  *)  R(l),  R(2),  R(3) 

C  RR=  SQRT  (R(l)**2  +  R(2)**2  +  R(3)**2) 

DO  351=1,3 
RH(I)=  R(I)  /  RR 
35  CONTINUE 

C  COMPUTE  REQUIRED  GONIOMETER  MIRROR  VECTOR  FROM  EQS  1-3 
DO  40 1=1,3 
N(I)=  (RH(I)  +  SH(I))  /  2 
40  CONTINUE 

NN=  SQRT(N(I)**2  +  N(2)**2  +  N(3)**2) 

DO  45 1=1,3 
NFa)=  N(I)  /  NN 
45  CONTINUE 
DL=  0 
PH=  0 

SD=  SIN  PL) 

CD=  COS  PL) 

SP-*  SIN  (PH) 

CP=  COS  (PH) 

C 

C  APPLY  INVERSE  OF  TRANSFORMATION  IN  EQS  8, 10 
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MH(1)=  NH(1)  *  CP  -  NH(3)  *  SP 

MH(2)=  NH(1)  *  SP  *  SD  +  NH(2)  *  CD  +  NH(3)  *  CP  *  SD 

MH(3)=  NH(1)  *  SP  *  CD  -  NH(2)  *  SD  +  NH(3)  *  CP  *  CD 


50 

C 

C 

11 

12 


C 

C 


C 

C 

C 

C 


C 


COMPUTE  GONIOMETER  ANGLES  TO  GIVEN  TARGET  COORDINATES 


M(l)--  MH(1) 

Mfl)=  MH(2) 

M(3)=  MH(3) 

MM=  SQRT  (M(l)**2  +  M(2)**2  +  M(3)**2) 
DO  50  ,3 

MH(I)=  M(I)  /  MM 
CONTINUE 


WRITE  (*,  108) 

READ  (10,  *1  R(3) 

FORMAT  ( K) 

FORMAT  (12) 

RETURN 

END 

**************«***»***«**a**************«********************* 
********* ***************************************************** 


SUBROUTINE  POS  (3E,DN,XAXIS,YAXIS,GP,PL,DL, ENTRY) 
INTEGER  XAXIS,YAXIS,GP 

INTEGER  DN 

INTEGERS  TIMO,  RN,  SE,  ENTRY, IPH.IDL 

INTEGERS  PL,  DL,  PM,  DM.  PS,  DS,  PD,  DD 

INTEGER *2  IBWRT,  IBRD 

PLt  DL  ARE  PHI,  DELTA  OF  LOCUS 

PM,  DM  ARE  PKI,  DELTA  OF  MANUAL  POSITION 

PS,  DS  ARE  PHI,  DELTA  OF  SCANNED  MAXIMUM 

PD,  DD  ARE  PHI,  DELTA  DIFFERENCES 

LOGICAL  BRANCH 

CHARACTERS  SS,  US,  G$ 

CHARACTERS  PAUSE 

CHARACTERS  BUFX.  BUFY 
CHARACTERS  CONV 
LOGICAL  FLAG,  FLAF 

LOGICAL  OP 

CHARACTER‘10  DRT,  FVAR 
CHARACTERS  1  LIST 

INTEGER  RCL 

FLAF =.  FALSE. 

ASK  USER  TO  MANUALLY  POSITION  GONIOMETER 
WRITE  (•,  3) 

IF  (ENTRY.EQ.l)  THEN 

PAUSE 

GOTO  20 

END  IF 

CALL  IBRD  (GP,  PAUSE,  6) 
LIST=,#7/*997/Y//77/’&,//,REDr//’@’ 

CALL  IBWRT  (GP,  LIST,  11) 
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C  SCAN  ACROSS  SENSOR  TO  IDENTIFY  POSITION  OF  MAX  IRRADIANCE 
20  CALL  SCAN  (XAXIS,  YAXIS,  GP,  PD,  DD,  ENTRY) 

C  SCAN  RETURNS  VALUES  PS,  DS  DISTANCES  TO  MAXIMUM 
FLAG =. FALSE. 

C  NOW  MOVE  TO  MAXIMUM 

GO  TO  50 
55  CONTINUE 

C  IPH=-10 

C  IDL=  10 

C  S$=’-7/CHAR(13) 

C  U$ = ’ + 7/CH  AR(  1 3) 

C  G$=’G’//CHAR(t3) 

C  CALL  IBCLR  (XAXIS) 

C  CALL  IBCLR  (YAXIS) 

C  CALL  MOVE1  (IPH,  IDL,  XAXIS,  YAXIS,  S$,  T$,  US,  V$,  GS) 

CALL  IBCLR  (XAXIS) 

CALL  IBCLR  (YAXIS) 

CALL  IBWRTF  (XAXIS,  ’DISPLAY’) 

CALL,  IBWRTF  (YAXIS,  ’DISPLAY’) 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBWAIT  (YAXIS,  TIMO) 

CALL  IBRD  (XAXIS,  BUFX,  9) 

CALL  IBRD  (YAXIS,  BUFY,  9) 

RN=  DN 

IF  (SE.EQ.l)  WRITE  (2,  5,  REC-RN)  DN,  BUFX,  BUFY 
IF  (SE.EQ.2)  WRITE  (3,  5,  REC=RN)  DN,  BUFX,  BUFY 
IF  (SE.EQ.3)  WRITE  (4,  5,  REC=RN)  DN,  BUFX,  BUFY 
IF  (SE.EQ.4)  WRITE  (8,  5,  REC=RN)  DN,  BUFX,  BUFY 
IF  (SE.EQ.5)  WRITE  (9,  5,  REC=RN)  DN,  BUFX,  BUFY 
C  READ  CHAR  STRING  INTO  INTEGER  FORMAT 
WRITE  (CONV,  7)  BUFX 
READ  (CONV,  8)  PS 
WRITE  (CONV,  7)  BUFY 
READ  (CONV,  8)  DS 
WRITE  (*,  1003)  PS,  DS 
FLAG=.TRUE. 

FLAF=,  FALSE. 

PD=  PL  -  PS 
DD=  DL  -  DS 

50  IF  (PD.GT.-100. AND.PD.LT.  100)  BRANCH  =  .TRUE. 

IF  (PD.GE,  100. OR. PD. LE. -100)  BRANCH  =  .FALSE. 

IF  (PD  .LT.  0)  GO  TO  100 
GOTO  110 

100  S$=’-’//CHARvi3) 

GO  TO  150 

110  SS=’  +  7/CHAR(ii 

150  IF  (DD  .LT.  0)  GO  T^  700 
GO  TO  210 

200  U$='-7/CHAR(13) 

GO  TO  250 

210  U$=’  +  ’//CHAR(13) 
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250  G$=’G7/CHAR(13) 

IF  (FLAF)  GOTO  51 
IF  (FLAG)  THEN 
WRITE  (*,  6) 

IF  (ENTRY.EQ.l)  THEN 

PAUSE 

GOTO  21 

END  IF 

CALL  IBRD  (GP,  PAUSE,  6) 

LIST=  *#7/*997/* : 7/’77/*&7/*REDI7/’<®  ’ 

CALL  IBWRT  (GP,  LIST,  11) 

21  WRITE  (*,  1002)  PL,  PS,  PD 
WRITE  (*,  1002)  DL,  DS,  DD 
FLAF =. TRUE. 

GOTO  50 
END  IF 

51  IF  (BRANCH)  GO  TO  9 

CALL  MOVE2  (PD,DD,XAXIS,YAXIS,SS,T$,U$,V$,G$) 

IF  (FLAG)  GOTO  10 
GOTO  55 

9  CALL  MOVE1  (PD,DD,XAXIS,YAXIS,S$,TS,U$,V$,G$) 

IF  (FLAG)  GOTO  10 

GOTO  55 

2  FORMAT  (ID 

3  FORMAT  (’0MANUALLY  POSITION  GONIOMETER  TO  DESCRIBED  TARGET  ’/) 

5  FORMAT  (I5,2A9) 

6  FORMAT  (*  WAITING  TO  RETURN  TO  HOME  POSITION  */) 

7  FORMAT  (A9) 

8  FORMAT  (19) 

1000  FORMAT  (’  \2A9)  . 

1002  FORMAT  (‘  LOCUS:  15,  *  SENSOR:  \  15,  *  DIFFERENCE:  \  15) 

1003  FORMAT  (’  X  POS  OF  MAX:  \  15,  ‘  Y  POS  OF  MAX:  ’,15) 

10  RETURN 
END 

SUBROUTINE  MOVE!  aX,IY,XAXIS,YAXlS,S$,TS,U$,VS,G$) 

C  ROUTINE  TAKES  RErURNED  STEP  COUNTS  AND  DIRECTS  CONTROLLER 
C  TO  MOVE  GONIOMETER  TO  SELECTED  TARGET 
C  THE  FOLLOWING  IS  GPIB  PROFESSIONAL  FORTRAN  HEADER 
integer*2  ibclr 
integer*2  ibrd 

integer*2  ibwrt,  TIMO,  END 
integer *2  BIN,  XEOS,  REOS,  LF,  S,  CMPL 
INTEGERS  TNONE,  T10us,T3tes,T100us,T300u5 
INTEGER*2  Tlms,T3ms,T10ms,T30ms,T100ms 
INTEGERS  T300ms,Tls,T3s,TIOs,T30s 
INTEGER *2  T100s,T300s,T1000s 
REAL  XSTEPS,  YSTEPS,  PH,  DL 

INTEGER  XAXIS,  YAXIS,  IOUT 
INTEGERS  DN.  A,  L,  K,  I 
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INTEGER*2  IXX,  IYY,  IX,  IY 


CHARACTERS 

CHARACTERS 

CHARACTERS 

CHARACTERS 

CHARACTERS 

CHARACTERS 

CHARACTERS 


N$ 

S$,  T$,  U$,  VS,  GS,  BAD 
X4$,  Y4$ 

X3$,  Y3$ 

X$,  YS,  Xl$,  Yl$ 

X2$,  Y2S 
RATE 


c  GPIB  Commands:  values 

DATA  TIMO/Z’4000’/,  END/Z’20007 
c  Iberr  error  messages:  values 
c  EOS  mode:  bit  values 

data  BIN/Z* 10007, XEOS/Z’8007,REOS/Z’4007 
c  Timeout  values 

data  TNONE/0/,T10us/l  /,T30us/2/,T100us/3/,T300us/4/ 
data  Tlms/5/,T3ms/6/,T10ms/7/,T30ms/8/,T100ms/9/ 
data  T300ms/l0/,Tls/ll/,T3s/12/,T10s/13/T3Qs/14/ 
data  T100s/15/,T300s/16/,T1000s/17/ 
data  S/08/,  LF/10/ 

RATE=  ’RATE.DAT’ 

A=  DN 

I  FORMAT  (’  OUTPUT  DATA  FILE  NAME:  ’) 

II  FORMAT  (A) 

16  FORMAT  (14) 

17  FORMAT  (13) 

18  FORMAT  (12) 

19  FORMAT  (II) 

C  RATE  IS  THE  FILE  WHICH  CONTAINS  A,R,S,F  STATEMENTS 
CALL  IBCLR  (XAXIS) 

CALL  IBWRTF  (XAXIS, RATE) 

CALL  IBCLR  (YAXIS) 

CALL  IBWRTF  (YAXIS.RATE) 

C 

IXX=!ABS(IX) 

IYY=IABS(IY) 

N$a=’N’ 

IF  (IXX.LT.IOOO  .AND.  IXX  .GE.  100)  WRITE(X$,17)IXX 
IF  (IXX.GE.1000)  WR1TE(X$,16)  IXX 
IF  (IYY.LT.  1000  .AND.  IYY  .GE.  100)  WRITE(Y$,17)  IYY 
IF  (IYY.GE.10QO)  WRITE(YS,16)  IYY 
IF  (IXX.LT.100.AND.1XX.GE.10)  WRITE  (XS,18)  IXX 
IF  aYY.LT.100.AND.IYY.GE.l0)  WRITE  (YS,18)  IYY 
IF  0XX.LT.  10)  WRITE  (X$,  19)  IXX 
IF  0YY.LT.  10)  WRITE  (Y$,  19)  IYY 
IF  (IXX.EQ.O)  X$=*00’ 

IF  OYY.EQ.O)  Y$=’00’ 

X1$=NS//’  7/XS(l:3)//CHAR(13) 

X2$=N$//’  7/XS(l:4)//CHAR(13) 

X3S-NS//'  7/XS(  1 :2)//CH AR(  1 3) 

X4$as NS//'  7/X$(l:  1)//CHAR(13) 
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Y1$=N$//’  7/Y$(l:3)//CHAR(13) 

Y2$=N$//’  7/Y$(l:4)//CHAR(13) 

Y3$=N$//’  7/Y$(l  :2)//CHAR(13) 

Y4$=N$//’  7/Y$(l :  1)//CHAR(13) 

IF(IXX.LT.  100.AND.IXX.GE.10.AND.IYY.LT.  100. AND.IYY.GE.  10)  GOTO  800 
IF  (IXX.LT.  100.AND.IXX.GE.  10. AND.IYY.GE.  1000)  GO  TO  900 
IF0XX.LT.  100.  AND.IXX.GE.  10.  AND.IYY.GE.  100.  AND.IYY.LT.  KXX))GOTO  12 
IF  (DOC.LT.  10.AND.IYY.LT.  10)  GO  TO  1300 
IF  0XX.LT.  10.AND.IYY.GE.10.AND.IYY.LT.  100)  GO  TO  1400 
IF  (IXX.lt.  10.AND.IYY.GE.100.AND.IYY.lt.  1000)  GOTO  1500 
IF  0XX.LT.  10. AND.IYY.GE.  1000)  GOTO  1600 
IF  0XX.GE.1O.AND.DCX.LT.  100.AND.IYY.LT.  10)  GOTO  1700 
800  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X3$,5) 

CALL  IBWRT(XAXIS,S$,2) 

CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBW AIT  (XAXIS,  TIMO) 

CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS ,  Y3$,5) 

CALL  IBWRT(Y  AXIS,U$,2) 

CALL  IBWRT(Y  AXIS  ,GS,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

900  CALL  1BCLR(XAXIS) 

CALL  IB WRT(X AXIS  ,X3$  ,5) 

CALL  IBWRT(XAXIS,S$,2) 

CALL  IBWRT(XAXIS,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBW  AIT  (XAXIS,  TIMO) 

CALL  IBCLR(YAXIS) 

CALL  IBWRT(YAX!S,Y2$,7) 

CALL  IBWRT(YAXIS,US,2) 

CALL  IBWRT(YAXIS,GS,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

12  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X3$,5) 

CALL  IBWRT(XAXIS,S$,2) 

CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBCLR(YAXIS) 

CALL  IBWRT(YAXIS,Y1S,6) 

CALL  1BWRT(YAXIS,US,2) 

CALL  IB WRT(Y AXIS  ,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

1300  CALL  1BCLR(XAX!S) 

CALL  IBWRT(XAXIS,X4S,4) 

CALL  IBWRT(XAX1S,SS,2) 
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CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS,  Y4$,4) 
CALL  IBWRT(YAXIS,U$,2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1400  CALL  IBCLR(XAXIS) 

CALL  BWRT(XAXIS,X4S,4) 
CALL  IBWRT(XAXIS,S$f2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IB  WATT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS,  Y3$,5) 
CALL  IBWRT(YAXIS,U$,2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1500  CALL  IBCLR(XAXIS) 

CALL  IBWRT(X  AXIS  ,X4$,4) 
CALL  !BWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(YAXIS,Y1S,6) 
CALL  IBWRT(Y AXIS,U$,2) 
CALL  IBWRT(Y AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1600  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X4S,4) 
CALL  IBWRT(XAXIS  ,S$,2) 
CALL  IBWRT(XAXIS,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS,  Y2S.7) 
CALL  IBWRT(YAXIS,US,2) 
CALL  IBWRT(YAX1S.GS,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1700  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X3$,5) 
CALL  IBWRT(XAX!S,S$,2) 
CALL  IBWRT(XAXIS,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
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CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS ,  Y4$,4) 

CALL  IBWRT(Y  AXIS»U$,2) 

CALL  IBWRT  (Y  AXIS ,  GS  ,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

300  CONTINUE 

RETURN 
END 

Q  ft************************************************** 

Q  *************************************************** 

SUBROUTINE  MOVE2  (IX,IY,XAXIS,YAXIS,S$,TS,U$,V$,G$) 

C  ROUTINE  TAKES  RETURNED  STEP  COUNTS  AND  DIRECTS  CONTROLLER 
C  TO  MOVE  GONIOMETER  TO  SELECTED  TARGET 

C  THE  FOLLOWING  IS  GPIB  PROFESSIONAL  FORTRAN  HEADER 

integer*2  ibclr 
integer*?  ibrd 
integer*2  ibwrt 

integer*2  BIN,  XEOS,  REOS,  TIMO,  END,  S,  LF 
INTEGERS  TNONE,  T10us,T30us,T100us,T300us 
INTEGERS  Tlms,T3im,T10ms,T30ms,T100ms 
INTEGERS  T300ms,Tls,T3s,T10s,T30s 
INTEGER *2  TIOOs,T300s,T1000s 
REAL  XSTEPS,  YSTEPS,  PH,  DL 

INTEGER  XAXIS,  YAXIS,  IOUT 
INTEGER*?  DN,  A,  L,  K,  J 
INTEGERS  IXX,  IYY,  IX,  IY 
CHARACTER*!  NS 

CHARACTER*2  S$,  TS,  US.  VS,  G$,  BAD 
CHARACTER*4  X4S,  Y4S 

CHARACTER’S  X3S,  Y3S 

CHARACTERS  XS,  Y$,  Xl$,  Y1S 
CHARACTER*?  X2S,  Y2S 

CHARACTERS  RATE 

c  GPIB  Commands:  values 

DATA  TIMO/Z’40007,  END/Z’20007 
c  Iberr  error  messages:  values 
c  EOS  mode:  bit  values 

data  BIN/Z* 10007,XEOS/Z’8007,REOS/Z*4007 
c  Timeout  values 

data  TNONE/O/.T  1  Ous/ 1  /,T30us/2/,T  1 0Ous/3/,T30Qus/4/ 
data  T1  ms/5/,T3ms/6/,Tl  0msG7,T30ms/8/,T  1 00tns/9/ 
data  T300ms/10/,Tls/l  l/,T3s/12/,T10s/13/T30s/I4/ 
dau  T100s/15/,T300s/16/,T!QOQs/l7/ 
dau  S/08/,  LF/10/ 

RATE=  'RATE.DAT' 

A=  DN 

11  FORMAT  (A) 

16  FORMAT  fl4) 

17  FORMAT  (U) 
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18  FORMAT  (12) 

19  FORMAT  ai) 

C  RATE  IS  THE  FILE  WHICH  CONTAINS  A,R,S,F  STATEMENTS 

CALL  IBCLR  (XAXIS) 

CALL  IBWRTF  (XAXIS, RATE) 

CALL  IBCLR  (YAXIS) 

CALL  IBWRTF  (YAXIS, RATE) 

C 

IXX=IABS(IX) 

IYY=IABS(IY) 

N$=’N* 

IF  (IXX.LT.  1000  .AND.  IXX  .GE.  100)  WRITE(X$,17)  DOC 

IF  (IXX.GE.  1000)  WRITE (X$,  16)  DOC 

IF  (IYY.LT.  1000  .AND.  IYY  .GE.  100)  WRITE(Y$,17)  IYY 

IF  (IYY.GE.1000)  WRTTE(Y$,16)  IYY 

IF  (IXX.LT.100.AND.DCX.GE.10)  WRITE  (X$,18)  DCX 

IF  0YY.LT.  100.  AND.DOC.GE.  10)  WRITE  (Y$,18)  IYY 

IF  (DCX.LT.  10)  WRITE  (X$,  19)  IXX 

IF  0YY.LT.10)  WRITE  (Y$,  19)  IYY 

IF  (IXX.EQ.0)  X$=W 

IF  0YY.EQ.O)  Y$=W 

X1$=N Sir  7/X$(l:3)//CHAR(13) 

X2$=N$//’  7/X$(l:4)//CHAR(13) 

X3$=N$//’  7/X$(l  :2)//CHAR(13) 

X4$=sN$//’  7/X$(  1 : 1  )//CH AR(  1 3) 

Y1S*NS//’  7/Y$(  1 :3)//CH AR(  1 3) 

Y2$=N$//‘  7/Y$(l  :4)//CHAR(13) 

Y3$«N$//’  7/Y$(  1 :2)//CH AR(  1 3) 

Y4$^^N$//,  7/Y$(l :  1)//CHAR(13) 

IF  0XX.LT.  1000.  AND.IXX.GE.  100. AND  JYY.GE.  1000)  GO  TO  500 

IF  OXX.GE.1000.AND.IYY.lt.  1000, AND.IYY.GE.  100)  GO  TO  600 

IF  0XX  .GE.  1000  .AND.  IYY  .GE.  1000)  GO  TO  700 

IF0XX.GE.  100.AND.IXX.LT.  1000.AND.IYY.LT.  100.AND.IYY.GE.  10)GOTO  10 

IF  (IXX.GE.  1000.AND.1YY.LT.  100. AND.IYY.GE.  10)  GOTO  1100 

IF  (IXX.GE. lOO.AND.IXX.LT.  1000.AND.IYY.LT.  10)  GOTO  1800 

IF  (IXX.GE.  lOOO.AND.IYY.LT.  10)  GO  TO  1900 

IF0XX.LT.  1000.  AND.DCX.GE.  100.  AND.IYY.LT.  1000.  AND.IYY.GE.  100) 

/  GOTO  20 

C  X  AXIS  MOVEMENT 

500  CALL  IBCLR(XAXIS) 

CALL  IBWRT(X  AXIS.X  1  $.6) 

CALL  !BWRT(XAXIS,S$.2) 

CALL  !BWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  1BCLR(YAXIS) 

CALL  1BWRT(YAX1S,Y2S,7) 

CALL  IBWRT(YAXIS,U$,2) 

CALL  IBWRT(Y  AXIS,GS,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
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600  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X2$,7) 
CALL  mWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  mWRT(YAXIS,YlS,6) 
CALL  IBWRT(Y  AXIS,U$,2) 
CALL  IBWRT(YAXIS,G$,2) 

2  GO  TO  300 

C  X  AXIS  MOVEMENT 

700  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X2$,7) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(Y  AXIS) 

CALL  IBWRT(Y  AXIS,Y2$,7) 
CALL  IBWRT(Y  AXIS, US, 2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
10  CALL  1BCLR(XAXIS) 

CALL  IBWRT(XAXIS,X1$,6) 
CALL  1BWRT(XAXIS,SS,2) 
CALL  IBWRT(X  AX  IS  ,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMC) 
CALL  IBCLR(Y  AXIS) 

CALL  IBWRT(YAXIS,Y3$,5) 
CALL  IBWRT(YAXIS,US,2) 
CALL  IBWRT(YAXIS,GS,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1100  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X2$,7) 
CALL  IBWRT(XAXIS,SS,2) 
CALL  IBWRT(XAXIS,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAX!S) 

CALL  IBWRT(YAXIS,Y3S,5) 
CALL  IBWRT(Y  AXIS, US, 2) 
CALL  IBWRT(Y  AXIS, GS, 2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1800  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X1$,6) 
CALL  IBWkT(XAXIS,SS,2) 
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CALL  IB WRT (X AXIS ,  G$ ,  2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBCLR(YAXIS) 

CALL  IBWRT(YAXIS,Y4$,4) 

CALL  IBWRT(Y AXIS,U$,2) 

CALL  IBWRT(Y AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1900  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X2$,7) 

CALL  IBWRT(XAXIS,S$,2) 

CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  1BCLR(YAXIS) 

CALL  IBWRT(YAXIS,Y4$,4) 

CALL  IBWRT(YAXIS,U$,2) 

CALL  IBWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
20  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X1$,6) 

CALL  IBWRT(XAXIS,S$,2) 

CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y AXIS, Yl$,6) 

CALL  IBWRT(YAXIS,U$,2) 

CALL  IBWRT(Y AXIS,GS,2) 

300  CONTINUE 
RETURN 
END 

c 

c 

SUBROUTINE  SCAN  (XAXIS,  YAXIS,  GP,  IX,  IY.  ENTRY) 

C  ROUTINE  STEPS  BEAM  THROUGH  SCANNING  OF  SENSOR  AND  DETERMINES 
C  POSITION  OF  MAXIMUM  IRRADIANCE 
INTEGER  XAXIS.  YAXIS,  GP 

INTEGER-2  XAR,  VAR.  MXAR,  MYAR 
INTEGER-2  1PH,  IDL,  IX,  IY.  DUM 
INTEGER-2  IBTMO,  1BCLR,  IBRD,  TIMO,  IBWRT 
INTEGER-2  T3MS.T3S.SPR 
INTEGER-2  ITYPE,  ENTRY 
REAL  ARRAY,  POWVAL 

CHARACTER-2  SS,  US,  GS 

CHARACTER-4  MESSAGE 

CHARACTER-6  POWER 

CHARACTER- 11  LIST 

DIMENSION  ARRAY  (1:5,  1:5) 
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DATA  T3MS/6/,TlS/l  1/.T3S/12/ 

C  PERFORM  DEVICE  CHECK 
CALL  IBRSP  (GP,  SPR) 

WRITE  (*,  *)  SPR 
IF  (SPR.NE.O)  GO  TO  90 

C  MOVE  BEAM  FROM  MANUALLY  SELECTED  POSITION  TO  BEGINNING  OF  SCAN 
IPH=0 
IDL=0 

SS=’-’//CHAR(13) 

US=’+7/CHAR(13) 

GS=’G7/CHAR(13) 

CALL  IBCLR  (XAXIS) 

CALL  IBCLR  (YAXIS) 

CALL  MOVE1  (IPH,  IDL,  XAXIS,  YAXIS,  S$,  T$,  U$,  V$,  GS) 

IF  (ENTRY.EQ.  1)  GOTO  31 
CALL  IBRD  (GP,  POWER,  6) 

WRITE  (*,  70)  POWER 

CALL  DECODE  (POWER,  ITYPE,  XAR,  YAR,  DN) 

MESSAGE  =  ’REDI’ 

IF  ((XAR.NE.  l).OR.(Y AR.NE.  1))  MESSAGE =’B  ADD* 

WRITE  (*,  *)  ’XAR  AND  YAR  ARE:’ 

WRITE  (*,  •)  XAR, YAR 

CALL  ENCODE  (XAR,  YAR,  ITYPE,  MESSAGE,  LIST) 

CALL  IBWRT  (GP,  LIST,  11) 

31  U$=’-7/CHAR(13) 

C  SCAN  THROUGH  25  STEPS  IN  5  INCREMENTS  OF  5  STEPS  EACH 
DO  10  J=5 
DO  20 1*1 

IF  (ENTRY.EQ.  1)  GOTO  32 
CALL  IBRD  (GP,  POWER.  6) 

CALL  DECODE  (POWER,  ITYPE,  XAR,  YAR,  DN) 

MESSAGE*  ’REDI' 

32  k=  m 

IF  (MOD(J,2).NE.Q)  THEN 
SS=’+’//CKAR(13) 

M=I+I 
END  IF 

IF  (MOD(J,2).EQ.O)  THEN 
SS=’-’//CHAR(l3) 

M=M 
END  IF 
N=J/5 
IPH=0 
IDL=0 

IF  (ENTRY.EQ.  1)  GOTO  33 
IF  (XAR.NE.M.OR.YAR.NE.N)  MESSAGE*  ’BADD’ 

33  CALL  MOVE1  (IPH.  IDL,  XAXIS,  YAXIS,  S$,  TS,  US,  V$,  GS) 

IF  (ENTRY.EQ.  1)  GOTO  20 

CALL  ENCODE  (XAR,  YAR,  tTYPE,  MESSAGE,  UST) 

CALL  IBWRT  (GP,  UST,  11) 

20  CONTINUE 
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IF  (J.EQ.5)  GOTO  II 

IPH=0 

IDL=0 

IF  (ENTRY.EQ.l)  GOTO  34 
CALL  IBRD  (GP,  POWER,  6) 

CALL  DECODE  (POWER,  ITYPE,  XAR,  YAR,  DN) 

MESSAGE®  ’REDF 

34  N=  N+l 

IF  (ENTRY.EQ.l)  GOTO  35 
IF  (XAR.NE.M.OR.YAR.NE.N)  MESSAGE®  ’BADD’ 

35  CALL  MOVE1  (IPH,  IDL,  XAXIS,  YAXIS,  S$,  TS,  US,  VS,  GS) 

IF  (ENTRY.EQ.  1)  GOTO  10 

CALL  ENCODE  (XAR,  YAR,  ITYPE,  MESSAGE,  US'D 
CALL  IBWRT  (GP,  UST,  11) 

10  CONTINUE 

11  IF  (ENTRY.EQ.l)  THEN 
MXAR®  1 

MYAR®  1 
GOTO  36 
END  IF 

CALL  IBRD  (GP.  POWER,  6) 

CALL  DECODE  (POWER,  ITYPE,  MXAR,  MYAR,  DN) 

IF  (XAR.LT.  1 .0R.XAR.GT.5.0R. YAR.LT.  1  .OR. YAR.GT.5)  MESSAGE® ’BADD’ 
IF  (ITYPE.GT.9.0R.ITYPE.LT.8)  MESSAGE®  ’BADD* 

IF  (ITYPE.EQ.9)  ENTRY® 3 

CALL  ENCODE  (MXAR,  MYAR,  ITYPE,  MESSAGE,  USD 
CALL  IBWRT  (GP,  LIST,  11) 

36  IX®  -(5*(1-MXAR)) 

IY®  (5M-MYAR)) 

RETURN 

70  FORMAT  (A  10) 

80  FORMAT  (110) 

90  WRITE  (*,•) 1  DEVICE  CHECK  ERROR' 

END 

SUBROUTINE  DECODE  (POWER,  ITYPE,  IDX,  INY,  DN) 

INTEGERS  ITYPE,  IDX,  INY,  DN 
CHARACTER* I  TYPE,  DX,  NY 

CHARACTERS  POWER,  CONV 

WRITE  (20.  1)  POWER 
TYPE®  POWER  (2:2) 

DX®  POWER  (4:4) 

NY®  POWER  (5:5) 

WRITE  (CONV,  2)  TYPE 
READ  (CONV,  3)  ITYPE 
WRITE  (*,  *)  ITYPE 
WRITE  (CONV,  2)  DX 
READ  (CONV,  3)  IDX 
WRITE  (*,  *)  IDX 
WRITE  (CONV,  2)  NY 
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READ  (C0NV,3)  INY 
WRITE  (*,  *)  INY 
DN=IDX*10 
DN=DN+INY 
WRITE  (*,  *)  DN 

1  FORMAT  (A6) 

2  FORMAT  (A) 

3  FORMAT  (II) 

RETURN 

END 

q  «*«*«************#**«*****«********************»******a* 

SUBROUTINE  ENCODE  (IPHI,  IDEL,  ITYPE,  MESSAGE,  LIST) 

INTEGERS  IPHI,  IDEL,  ITYPE 

CHARACTER*  1  PHI,  DEL,  TYPE,  CONV 

CHARACTER*4  MESSAGE 

CHARACTER*1 1  LIST 

WRITE  (CONV,  1)  IPHI 

READ  (CONV ,2)  PHI 

WRITE  (CONV,  1)  IDEL 

READ  (CONV,  2)  DEL 

WRITE  (CONV,  1)  ITYPE 

READ  (CONV,  2)  TYPE 

LIST-  ,#7/PHI//DEL//*:7/TYPE//'&7/MESSAGE//,@’ 

WRITE  (*,  2)  UST 

1  FORMAT  (11) 

2  FORMAT  (A) 

RETURN 

END 
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C2.  Pointing  Routine 


C 

C 

C 

C 

C 

C 

C 

C 

C 


FROM  PROGRAM  HITRGT.FOR  SEPT  15,  1987 
LINK  WITH  PFIB; 

PROGRAM  TO  READ  TARGET  POSITION  DATA  FILES: 

OAR.DAT,  UVH.DAT,  UVL.DAT,  IRH.DAT,  IRLDAT, 

FROM  DATA  RECEIVED  FROM  LACE  MVAX,  AND  TO  THEN 
EXECUTE  POSITIONING  OF  LASER  BEAM  ON  GIVEN  TARGET 
PROGRAM  ALSO  TAKES  POWER  MEASUREMENT,  CALIBRATES  AND 
RETURNS  POWER  TO  LACE  MVAX. 

GPIB  PROFESSIONAL  FORTRAN  HEADER 
COMMON  /IBGLOB/  IBSTA,  IBERR,  IBCNT 
GPIB  FUNCTIONS 

INTEGERS  IBSTA,  IBERR,  IBCNT 
INTEGER*2  IBCLR,  IBCMD,  IBEOS,  EBFIND,  IBONL 
INTEGER*2  IBRD,  IBRDF,  IBRSP,  IBSRE,  IBTMO 
INTEGERS  IBWAJT,  IBWRT,  IBWRTF 
GPIB  COMMANDS 

INTEGERS  UNL,  UNT,  GTL,  SDC,  PPC,  GET,  TCT,  SPR 
INTEGER*2  LLO,  DCL,  PPU,  SPE,  SPD,  PPE,  PPD 
INTEGER*2  ERR,  TIMO,  END,  SRQI,  RQS,  CMPL,  LOK 
INTEGER*2  REM.  CIC,  ATN,  TACS,  LACS,  DTAS,  DC  AS 
TIMEOUT  VALUES 

INTEGER*2  TNONE,  T10US,  T30US,  T100US,  T300US 
INTEGER*2  TIMS,  T3MS,  T10MS,  T30MS,  T100MS 
INTEGER*2  T300MS,  T1S,  US,  T10S,  T30S 
IN1EGER*2  T100S,  T300S,  T10C0S 
MAIN  PROGRAM  VARIABLE  DECLARATIONS 
INTEGER  XAXIS,  YAXIS,  ADC,  GP 

INTEGER*?  PH,  DL,  IX,  IY,  DN,  DATA 
INTEGERS  ITYPE,  VAL,  LASER,  LS 
REAL  POUT,  BEAM,  PRF,  PL,  RATIO 


BUF 

TYPE,  CONV 
CDN,  DONV 
BONV, CPOUT 
FNAME 
LIST 


C 

C 


CHARACTER*6 
CHARACTER*  1 
CHARACTER*2 
CHARACTER*? 

CHARACTER*  10 
CHARACTER*  14 
LOGICAL  OP,  EX 
GPIB  COMMANDS:  VALUES 

DATA  UNL/C3/,  UNT'95/,  GTL/01/,  SDC/04/,  PPC/05/ 
DATA  GET/08/,  TCT/09/,  LLO/17/,  DCL/20/,  PPU/21/ 

DATA  SPE/24/,  SPD/25/,  PPE/96/,  PPD/112/ 

GPIB  STATUS  BIT  VECTOR:  VALUES 
DATA  TIMO/Z’ 40007,  END/Z’20007,  SRQI/Z’ 10007 
DATA  CIC/Z’20’/,  RQS/Z’8007.  CMPL/Z7007,  LOK/Z’807 
DATA  REM/Z’407,  ATN/Z’107,  TACS/Z’87,  LACS/Z’47 
DATA  DTAS/Z’27,  DCAS/ZT/ 

EOS  MODE:  BIT  VALUES 

DATA  BIN/Z’  10007, XEOS/ZW/,REOS/Z’4007 

TIMEOUT  VALUES 
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DATA  TNONE/0/,T10us/l/,T30us/2/,T100us/3/,1300us/4/ 

DATA  T1  ms/5/,T3ms/6/,T  1 0ms/7/,T30ms/8/,T  1 00ms/9/ 

DATA  T300ms/10/,Tls/ll/,T3s/12/T10s/13/,T30s/14/ 

DATA  T100s/15/,T300s/16/,T1000s/17/ 

C  DEVICE  ADDRESS  LOCATION 

XAXIS  =  IBFIND  fX  AXIS’) 

YAXIS=  IBFIND  (TAXIS’) 

GP=  IBFIND  (’GPIB2’) 

ADC=  IBFIND  (’ADC’) 

C  PERFORM  DEVICE  CONNECTION  CHECKS 

CALL  IBRSP  (XAXIS,  SPR) 

WRITE  (*,  *)  SPR 
IF  (SPR.NE.48)  GOTO  12 
CALL  IBRSP  (YAXIS,  SPR) 

WRITE  (*,  *)  SPR 
IF  (SPR.NE.48)  GOTO  12 
CALL  IBRSP  (GP,  SPR) 

WRITE  (*  *)  SPR 
CALL  IBRSP  (ADC,  SPR) 

WRITE  (*,  *)  SPR 

C  AN  IBWRT  WILL  HANG  ON  THE  BUS  UNTIL  MVAX  READ  ADDRESSES 
C  IBMPC  TO  WRITE  TO  BUS 

C  AN  IBRD  WILL  HANG  ON  THE  BUS  UNTIL  MVAX  WRITES  DATA  ON  BUS 
C  ASK  USER  WHETHER  TO  EXPECT  DATA  FROM  BUS  OR  KEYBOARD 
WRITE  (*,*)  ’EXPECT  DATA  FROM  KEYBOARD  (1),  OR  IEEE  BUS  (2)’ 
READ  (*,  *)  DATA 
1  CONTINUE 

CALL  SEARCH  (PH,  DL,  DN,  ITYPE,  DATA,  GP) 

C  EXECUTE  MOVEMENT 

CALL  MOVE  (PH,  DL,  XAXIS,  YAXIS) 

C  TAKE  POWER  MEASUREMENT  READING  FROM  ADC 

CALL  POWER  (ADC,  ITYPE,  POUT) 

IF(DATA.EQ.l)  GOTO  1 

C  ENCODE  POUT  AS  DATA  WORD  TO  SEND  BACK  TO  MVAX 
WRITE  (DONV,  20)  DN 
READ  (DONV,  21)  CDN 
WRITE  (CONV,  22)  ITYPE 
READ  (CONV,  23)  TYPE 
WRITE  (BONV,  24)  POUT 
READ  (BONV,  25)  CPOUT 
UST=  *#’//CDN//’ :  7/TYPE//’&7/CPOUT//’@  ’ 

C  SEND  POWER  MEASUREMENT  TO  MVAX 
WRITE  (•,  26)  LIST 
CALL  IBWRT  (GP,  LIST,  14) 

C  LOOP  BACK,  GET  ANOTHER  SENSOR  ID  NUMBER 
GOTO  1 

20  FORMAT  (12.2) 

21  FORMAT  (A2) 

22  FORMAT  (H) 

23  FORMAT  (A  1) 

24  FORMAT  (F7.3) 
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25  FORMAT  (A7) 

26  FORMAT  ( A) 

11  FORMAT  (A5) 

12  WRITE  (*  *)  ’DEVICE  CHECK  FAILURE,  POWER  DOWN  CONTROLLER’ 
END 

C  ********************************************************** 

SUBROUTINE  SEARCH  (PH,  DL,  DN,  ITYPE,  K,  GP) 

C  ROUTINE  DECODES  DATA  READ  FROM  MVAX,  OPENS  DATA  FILE, 

C  RETURNS  WITH  DN,  PH,  DL  VALUES 

INTEGER*2  PH,  DL,  DN,  DM 
INTEGERS  ITYPE,  K,  IBRD,  IBWRT 
INTEGER  RL,  CODE,  GP 

CHARACTER*?  FILNM 

CHARACTER*6  BUF 
CHARACTER*4  DONV 

CHARACTER*  1  CONV 

CHARACTER*  1  BDN,  BDM 

CHARACTER*!  BTYPE 

RL=  30 

103  IF  (K.EQ.l)  THEN 
WRITE  (*,  99) 

WRITE  (*,  101) 

READ  (*,  *)  ITYPE,  DN 
GO  TO  102 
END  IF 

IF  (K.EQ.2)  THEN 
CALL  IBRD  (GP,  BUF,  6) 

BTYPE  =  BUF(2:2) 

BDN*  BUF(4:4) 

BDM*  BUF(5:5) 

WRITE  (CONV,  4)  BTYPE 
READ  (CONV,  6)  ITYPE 
WRITE  (CONV,  4)  BDN 
READ  (CONV,  6)  DN 
WRITE  (CONV,  4)  BDM 
READ  (CONV,  6)  DM 
DN*DN*10 
DN*DN  +  DM 
END  IF 

102  CODE*  85 

IF  (ITYPE.EQ.4.0R.ITYPE.EQJ)  CODE*  40 
IF  (ITYPE.EQ.l)  FILNM* ’OAR. DAT* 

IF  (ITYPE.EQ.2)  FILNM  =  ’UVH.DAT’ 

IF  (ITYPE.EQ.3)  FILNM =’UVL.DAT’ 

IF  (ITYPE.EQ.4)  FILNM  =  ’IRH.DAT’ 

IF  (ITYPE.EQ.5)  FILNM  =  ’IRL.DAT’ 

IF  (ITYPE.LT. l.OR.ITYPE.GT.5)  GO  TO  3 
WRITE  (*,  •)  FILNM 

OPEN  (2,  FILE=FILNM,  ACCESS  *  ’DIRECT’,  RECL=RL, 

I  FORM*  ’FORMATTED’,  STATUS=’OLD’) 

C  SEARCH  FOR  PH,  DL 
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10 

1 

3 

4 

5 

6 

7 

8 
9 

99 

100 

/ 

101 

c 

c 

c 


c 


c 


c 


c 


READ  (2,  REC=DN,  FMT=8)  J,  PH,  DL 

WRITE  (*,  *)  J,  PH,  DL 

CONTINUE 

RETURN 

WRITE  (*,  9) 

GOTO  103 
FORMAT  (A) 

FORMAT  (A3) 

FORMAT  (II) 

FORMAT  (A4) 

FORMAT  (15,  219) 

FORMAT  (’  SENSOR  TYPE  NOT  PROPERLY  IDENTIFIED  ’) 

FORMAT  (’  SENSOR  TYPES:  (1)AR,  (2)UVH,  (3)UVL,  (4)IRH,  (5)IRL’) 
FORMAT  C  ENTER  (1)  FOR  USER  INPUT, 

(2)  FOR  IEEE  DATA  XFER:  *) 
FORMATf  ENTER  SENSOR  TYPE,  ID#:  ’) 

END 

********************************************************* 


SUBROUTINE  MOVE  (PH,  DL,  XAXIS,  YAXIS) 

ROUTINE  TAKES  RETURNED  STEP  COUNTS  AND  DIRECTS  CONTROLLER 
TO  MOVE  GONIOMETER  TO  SELECTED  TARGET 
INTEGERS  PH,  DL,  IXX,  IYY,  IX,  IY 
INTEGER*2  XAXIS,  YAXIS,  T1MO 


CHARACTER*  1 
CHARACTER*2 
CHARACTER*4 
CHARACTER*5 
CHARACTER*6 
CHARACTER*? 
CHARACTER'S 
CHARACTER*9 


NS 

S$,  US,  GS 

X4S,  Y4S 

X3S,  Y3S 

X$,  Y$,  X1S,  Y1S 

X2S,  Y2S,  BUFX,  BUFY 

RATE 


CV 


READ  CURRENT  POSITION  FROM  DISPLAY 
CALL  IBCLR  (XAXIS) 

CALL  IBCLR  (YAXIS) 

CALL  IBWRTF  (XAXIS,  ’DISPLAY’) 

CALI..  IBWRTF  (YAXIS,  ’DISPLAY*) 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBWAIT  (YAXIS,  TIMO) 

CALL  IBRD  (XAXIS,  BUFX,  9) 

CALL  IBRD  (YAXIS,  BUFY,  9) 

WRITE  (*,  2000)  BUFX,  BUFY 

READ  CHARACTER  STRING  INTO  INTEGER  FORMAT 


WRITE  (CV,  70)  BUFX 
READ  (CV,  80)  IX 
WRITE  (CV,  70)  BUFY 
READ  (CV,  80)  IY 

SUBTRACT  PRESENT  VALUES  FROM  DESTINATION  POSITION 


PH-  PH  -IX 
DL=  DL  -  IY 
RATE=  ’RATE.DAT’ 

RATE  IS  THE  FILE  WHICH  CONTAINS  A,R,S,F  STATEMENTS 
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CALL  IBCLR  (XAXIS) 

CALL  IBWRTF  (XAXIS, RATE) 

CALL  IBCLR  (YAXIS) 

CALL  IBWRTF  (YAXIS, RATE) 

C  WRITE  STEP  COUNTS  AS  STRINGS 
IXX=IABS(PH) 

IYY=IABS(DL) 

N$=’N’ 

IF  GXX.LT.  1000  .AND.  DOC  .GE.  100)  WRITERS,  17)  DOC 

IF  (IXX.GE.1000)  WRTTE(X$,16)  IXX 

IF  (IYY.LT.  1000  .AND.  IYY  .GE.  100)  WRITE(YS,17)  IYY 

IF  (IYY.GE.1000)  WRTTE(YS,16)  IYY 

IF  (IXX.LT.  100. AND. IXX. GE.  10)  WRITE  (X$,18)  IXX 

IF  GYY.LT.100.AND.IYY.GE.10)  WRITE  (Y$,18)  IYY 

IF  GXX.LT.  10)  WRITE  (XS,  19)  DOC 

IF  (IYY.LT.  10)  WRITE  (Y$,  19)  IYY 

IF  (IXX.EQ.0)  XS= W 

IF  GYY.EQ.0)  YS=W 

X1$=NS//’  ’//X$(l  :3)//CHAR(13) 

X2$=N$//’  7/X$(l  :4)//CHAR(13) 

X3$=*N$//’  7/XS(  1 :2)//CH AR(  1 3) 

X4$=N$/'f  7/XS(  1 : 1  )//CH AR(  1 3) 

Y1$=N$//*  7/Y$(l  :3)//CHAR(13) 

Y2$=N$//’  7/YS(l  :4)//CHAR(13) 

Y3$=N$//*  7/Y$(l  :2)//CHAR(13) 

Y4S=N$//*  7/Y$(  1 : 1  )//CH AR(1 3) 

IF  (PH  XT.  0)  GO  TO  100 
GOTO  110 

100  S$=’-7/CHAR(13) 

GOTO  150 

110  SS=’+7/CHAR(13) 

150  IF  G)L  XT.  0)  GO  TO  200 
GO  TO  210 

200  U$=’-7/CHAR(13) 

GO  TO  250 

210  US=’+7/CHAR(13) 

250  G$=’G7/CHAR(13) 

IFGXX.LT.  100.  AND.IXX.GE.  10. AND.IYY.LT.  100.  ANDJYY.GE.  10)GOTO  800 

IF  GXX.LT.  100.  AND.IXX.GE.  10. AND.IYY.GE.  1000)  GO  TO  900 

IFGXX.LT.  100.  AND.IXX.GE.  10.  AND.IYY.GE.  100.  AND.IYY.LT.  1000)GOTO  12 

IF  GXX.LT.  10.AND.1YY.LT.  10)  GOTO  1300 

IF  GXX.LT.  10.AND.IYY.GE.10.AND.IYY.LT.  100)  GOTO  MOO 

IF  GXX.LT.  1O.AND.IYY.GE.I00.AND.IYY.LT.  1000)  GOTO  1500 

IF  GXX.LT.  10. AND.IYY.GE.  1000)  GO  TO  1600 

IF  GXX.GE.10.AND.IXX.LT.  100. AND.IYY.LT.10)  GO  TO  1700 

IF  GXX.LT.  1000. AND.IXX.GE.  100. AND.IYY.GE.  1000)  GO  TO  500 

IF  GXX.GE.  1000. AND.IYY.LT.  1000. AND.IYY.GE.  100)  GOTO  600 

IF  GXX  .GE.  1000  .AND.  IYY  .GE.  1000)  GO  TO  700 

IFGXX.GE.  100.AND.IXXXT.  1000.AND.IYYXT.  100.  AND.IYY.GE.  10)GOTO  10 

IF  GXX.GE.  1000.AND.IYY.LT.  100.AND.IYY.GE.  10)  GOTO  1100 

IF  GXX.GE.  100. AND.IXX.LT.  1000. AND.IYY.LT.  10)  GO  TO  1800 
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IF  (KX.GE.1000.AND.IYY.lt.  10)  GO  TO  1900 
IF(IXX.LT.  1000.AND.IXX.GE.  100.AND.IYY.LT.  1000.  AND.IYY.GE.  100) 
/  GOTO  20 

800  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X3$,5) 

CALL  BWRT(XAXIS,SS>2) 

CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBCLR(Y AXIS) 

CALL  IBWRT(Y  AXISfY3$,5) 

CALL  IBWRT(YAXIS,U$,2) 

CALL  IBWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

900  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X3$,5) 

CALL  IBWRT(XAXIS,S$,2) 

CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y AXIS,  Y2$,7) 

CALL  IBWRT(Y AXIS,U$,2) 

CALL  IBWRT(Y AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

12  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X3$,5) 

CALL  IBWRT(XAXIS,SS,2) 

CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y AXIS,Y  1  $,6) 

CALL  IBWRT(Y AXIS, US, 2) 

CALL  IBWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

1300  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X4$,4) 

CALL  IBWRT(XAXIS,S$,2) 

CALL  IBWRT(XAX1S,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 

CALL  IBCLR(Y  AXIS) 

CALL  IBWRT(Y AX1S,Y4$,4) 

CALL  IBWRT(YAXIS,U$,2) 

CALL  IBWRT(YAX1S,GS4) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
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1400  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X4$,4) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IB  WAIT  (X AXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(YAXIS,Y3$,5) 
CALL  IBWRT(Y  AXIS,U$,2) 
CALL  IBWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

1500  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXISfX4$,4) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS,Y1$,6) 
CALL  IBWRT(Y AXIS,U$,2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

1600  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X4$,4) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  1BWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  1BWRT(Y  AXIS,  Y2$,7) 
CALL  1BWRT(YAXIS,U$,2) 
CALL  IBWRT(Y AX1S,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

1700  CALL  1BCLR(XAXIS) 

CALL  IBWRT(XAXIS,X3$,5) 
CALL  IBWRT(XAXIS,SS,2) 
CALL  IBWRT<XAXI$,GS42) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y AXIS, Y4S.4 ) 
CALL  1BWRT(YAXIS,1)$,2) 
CALL  IBWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

500  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,Xl$,6) 
CALL  IBWRT(XAXIS,S$,2) 
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CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS,Y2$,7) 
CALL  IBWRT(YAXIS,U$,2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

600  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X2$,7) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(Y  AXIS) 

CALL  IBWRT(Y  AXIS,Y1$,6) 
CALL  IBWRT(Y AXIS,U$,2) 
CALL  IBWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

700  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X2$,7) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(Y  AXIS) 

CALL  IBWRT(YAXIS,Y2$,7) 
CALL  IBWRT(Y AXIS,U$,2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

10  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X1$,6) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS  ,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y AXIS,Y3$,5) 
CALL  IBWRT(Y AXIS,US,2) 
CALL  IBWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 

1100  CALL  IBCLR(XAX1S) 

CALL  IBWRT(XAX1S,X2$,7) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  1BWRT(XAX1S,GS,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
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CALL  mCLR(YAXIS) 

CALL  IBWRT(Y AXIS,Y3$,5) 
CALL  IBWRT(YAXIS,U$,2) 
CALL  mWRT(YAXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1800  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X1$,6) 
CALL  IBWRT(XAXIS,S$,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  DBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y  AXIS,Y4$,4) 
CALL  IBWRT(Y  AXIS,U$,2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
1900  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X2$,7) 
CALL  IBWRT(XAXIS,SS,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS,  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(YAXIS,Y4S,4) 
CALL  IBWRT(Y  AXIS, US, 2) 
CALL  IBWRT(Y  AXIS,G$,2) 

GO  TO  300 

C  X  AXIS  MOVEMENT 
20  CALL  IBCLR(XAXIS) 

CALL  IBWRT(XAXIS,X1S,6) 
CALL  IBWRT(XAXIS,SS,2) 
CALL  IBWRT(XAXIS,G$,2) 

C  Y  AXIS  MOVEMENT 

CALL  IBWAIT  (XAXIS.  TIMO) 
CALL  IBCLR(YAXIS) 

CALL  IBWRT(Y AXIS, Y  1S,6) 
CALL  IBWRT(Y  AXIS, US, 2) 
CALL  IBWRT(YAX1S,GS,2) 

300  CONTINUE 
11  FORMAT  (A) 

16  FORMAT  (M) 

17  FORMAT  (13) 

18  FORMAT  (12) 

19  FORMAT  (II) 

70  FORMAT  (A9) 

80  FORMAT  09) 

2000  FORMAT  (2A9) 

RETURN 

END 
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c 

c 

c 

c 

23 

C 

C 


21 

C 

22 


C 


C 

C 

C 

C 

C 


ft************:*************#**#**************************** 


SUBROUTINE  POWER  (ADC,  ITYPE,  POUT) 
INTEGERS  PIN,  ADC,  ITYPE,  LASER 


INTEGER  NUM 

REAL  RATIO,  SCALE,  ND 

REAL  BEAM,  PRF,  PL,  SENSAR,  AREA,  POUT 

CHARACTERS  BAD 

CHARACTERS!)  FILNM 

LOGICAL  EX,  OP 

CALL  IBCLR  (ADC) 

CALL  IBWRT  (ADC,  ’H8AJ\  4) 

CALL  IBRD  (ADC,  BAD,  2) 

CALL  CONVERT  (BAD,  PIN) 

WRITE  (*,  *)  PIN 


WRITE  (*,  5) 

READ  (*,  6)  SCALE 
SCALE  =  1.0 
WRITE  (*,  7) 

READ  (*,  8)  ND 
ND  =  1 
CONTINUE 
WRITE  (*,  10) 

READ  (*,  *)  LASER 
LASER  a  1 

OPEN  (9,  FILE=*LASERS.DAT\  STATUS =*OLD*) 
DO  21 1*1,5 

READ  (9,  *)  NUM,  BEAM,  PRF,  PL,  RATIO 
IF  (1.EQ.5.AND.NUM.NE.LASER)  GO  TO  23 
IF  (NUM.EQ.LASER)  GO  TO  22 
CONTINUE 

WHAT  IS  SENSOR  AREA? 

CLOSE  (9) 

IF  (ITYPE.EQ.l)  SENSAR*  1.0 
IF  (ITYPE.EQ.2)  SENSAR*  1.0 
IF  (ITYPE.EQ.3)  SENSAR*  1.0 
IF  (ITYPE.EQ.4)  SENSAR*  1.0 
IF  (ITYPE.EQ.5)  SENSAR*  1.0 
COMPUTE  FACTORS 
PIN*  PIN  +  1024 
SCALE*  SCALE  / 1024 


ND*  10  ••  (-ND) 

PRF*  1  /  PRF 
PL*  1/  PL 

AREA*  (SENSAR  /  BEAM) 

POUT*  PIN  •  SCALE  •  ND  •  PRF  •  PL  •  RATIO  •  AREA 
WRITE  (•,  •)  POUT 

FORMAT  (*  POWER  MEIER  FULL  SCALE  VALUE:  ') 

FORMAT  (F5.2) 

FORMAT  (’  NEUTRAL  DENSITY  VALUE: ') 

FORMAT  (F5.2) 

FORMAT  (’  SELECT  LASER:  (1)AR,  (2)HN,  (3)DF,  (4)EX,  (5)YAG  *) 
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RETURN 

END 

q  0t0m*^*****»i**********-''************************************ 

SUBROUTINE  CONVERT  (BAD.IOUT) 

C  THIS  ROUTINE  CONVERTS  THE  ADC  READING  TO  AN  INTEGER 
C  ROUTINE  IS  CALLED  UP  BY  SUBROUTINE  POWER 

CHARACTER*2  BAD 
CHARACTER*  1  BADA.BADB 
INTEGERS  IA,IB,IC,ID,IE,IGF,IOUT 
BADA=BAD(:1) 

BADB=BAD(2:) 

IA = ICH  AR(B  AD  A) 

IB = ICH  AR(BADB) 

IF  (IA  .LT.  4)  GO  TO  200 
100  IC  =  255  -  IA 
ID  »  IC*256 
IE  =  256  -  IB 
IGF  *  ID  +  IE 
IOUT  «  -1  *  IGF 
GO  TO  300 
200  IC  =  IA  *  256 
IOUT  =  IC  +  IB 
300  RETURN 
END 
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